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1.  Introduction 


Ceramics  have  several  properties  that  make  them  appealing  materials  for  a  wide  range 
of  practical  applications.  Since  their  strength  is  limited  by  brittle  failure,  the  fracture  of 
ceramics  has  been  the  subject  of  extensive  research.  However,  despite  recent  progress  in  the 
area,  a  number  of  issues  remain  unresolved.  For  example,  while  the  measured  toughness 
of  single  crystals  is  comparable  to  their  surface  energy,  that  of  polycrystalline  materials  is 
found  to  be  up  to  an  order  of  magnitude  greater.  In  addition,  the  apparent  toughness  of  the 
material  is  observed  to  increase  with  crack  length.  Three  mechanisms  have  been  proposed 
that  may  explain  this  behaviour:  crack  tip  shielding  by  microcracking;  crack  deflection; 
and  crack  trapping  combined  with  crack  face  bridging  by  tough  particles. 

The  possibility  that  microcracking  may  improve  toughness  has  recently  been  the  sub¬ 
ject  of  intense  investigation.  Extensive  microcracking  is  in  fact  observed  around  the  tips  of 
cracks  in  brittle  polycrystalline  materials  (Hoagland  et  al,  1973;  Claussen,  1976;  Mecholsky 
et  al,  1976;  Wu  et  al,  1978).  The  microcracks  increase  the  compliance  of  the  solid,  and 
so  shield  the  crack  tip  from  remote  stresses.  It  has  been  estimated  that  for  a  stationary 
mode  I  crack,  the  near-tip  stress  intensity  factors  are  reduced  by  around  30%  due  to  mi¬ 
crocracking.  (Green  et  al,  1973;  Hoagland  et  al,  1980;  Evans  and  Faber,  1981;  Evans,  1984; 
Evans  and  Fu,  1985;  Kachanov,  1986;  Hutchinson,  1987;  Rodin,  1987;  Ortiz,  1987,  1988; 
Ortiz  and  Giannakopoulos,  1989).  The  effect  is  slightly  greater  for  a  propagating  crack, 
due  to  the  influence  of  the  microcracked  wake.  (Charalambides  and  McMeeking,  1987; 
Ortiz  and  Giannakopoulos,  1989).  However,  for  a  crack  in  mixed  mode  or  mode  II,  the 
shielding  effect  is  almost  imperceptible  (Ortiz  and  Giannakopoulos,  1989).  Furthermore, 
although  the  microcracks  do  shield  the  crack  tip,  they  also  weaken  the  material  ahead  of 
the  crack.  Ortiz  (1988)  has  estimated  that  the  reduction  in  strength  due  to  microcracking 
almost  exactly  offsets  any  gain  in  touhness  due  to  crack  tip  shielding.  Microcracking  does 
not  therefore  appear  to  be  an  effective  toughening  mechanism. 

Crack  deflection  is  a  second  plausible  mechanism  for  toughening  brittle  materials. 
Since  the  fracture  of  ceramics  is  predominantly  intergranular,  cracks  tend  to  follow  a  tortu¬ 
ous  path  through  the  material.  Second  phase  particles  and  inter-granular  residual  stresses 
also  tend  to  deflect  the  path  of  a  crack.  This  increases  the  fracture  surface  area,  and 
consequently  raises  the  apparent  toughness  of  the  material.  A  rigorous  three-dimensional 
analysis  of  crack  deflection  has  not  been  attempted  to  date,  but  preliminary  estimates  sug¬ 
gest  that  crack  deflection  may  incu-us^.  toughness  by  around  10-20%  (Faber  and  Evans, 
1983,  Suresh.  1985). 
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Crack  trapping  and  crack  face  bridging  are  by  far  the  most  effective  toughening  mech¬ 
anisms.  Several  experiments  have  demonstrated  that  if  small  quantities  of  ductile  particles 
are  added  to  a  brittle  matrix  material,  its  toughness  may  be  increased  by  up  to  60  times 
(Mendelson  and  Fine,  1974;  Krstic  et  al,  1981).  It  has  long  been  recognised  that  crack 
bridging  is  responsible  for  improving  the  toughness  of  fibrc-reinforced  ceramics  and  ce¬ 
ment.  Bridging  particles  have  also  been  observed  in  the  wake  of  cracks  in  monolithic 
ceramics  (Swanson  et  al,  1987).  Fig.  1  illustrates  the  mechanism  for  forming  bridging 
particles  in  the  wake  of  a  crack.  As  a  semi-infinite  crack  propagates  through  a  material 
containing  a  distribution  of  particles,  parts  of  the  crack  front  reach  the  regions  of  higher 
toughness  and  consequently  arrest.  As  the  load  is  increased,  the  remainder  of  the  crack 
continues  to  propagate,  and  bows  out  between  the  pinning  particles.  This  process  of  crack 
trapping  is  an  important  source  of  toughness  in  itself  (Lange,  1970).  If  the  particles  are 
sufficiently  tough,  the  crack  front  eventually  bypasses  the  obstacles,  so  that  they  are  left 
behind  in  the  wake  of  the  crack. 

While  the  origin  of  bridging  particles  in  composite  materials  is  obvious,  it  is  not  imme¬ 
diately  apparent  why  bridging  particles  should  be  observed  in  monolithic  ceramics.  Evans 
et  al  (1977)  have  suggested  that  the  bridging  particles  may  be  attributed  to  the  action 
of  inter-granular  residual  stresses.  Substantial  residual  stresses  may  be  set  up  in  poly¬ 
crystalline  ceramics  during  manufacture,  due  to  the  thermal  and  elastic  anisotropy  of  the 
grains:  stresses  of  the  order  70-150  Mpa  have  been  measured  experimentally  (Blendel  and 
Coble,  1982)  and  stresses  of  a  similar  magnitude  have  been  predicted  theoretically  (Ortiz 
and  Molinari,  1988,  Ortiz  and  Suresh,  1990).  According  to  Evans  et  al  (1977),  grains  in  a 
state  of  residual  compression  are  less  likely  to  be  fractured  by  a  growing  crack  and  may  be 
left  behind  as  bridging  particles.  Evans  has  pointed  out  that  compressive  grains  may  be 
thought  of  as  regions  of  increased  toughness,  and  has  suggested  an  expression  relating  the 
magnitude  of  the  compressive  stress  and  the  grain  size  to  the  equivalent  toughness  of  the 
particles. 

Several  criteria  must  be  satisfied  in  order  to  maximise  the  toughening  due  to  crack 
bridging.  The  elastic  modulus  of  the  particles  must  not  exceed  that  of  the  matrix  material, 
otherwise  the  crack  tends  to  be  deflected  away  from  the  particles.  The  toughening  is  also 
sensitive  to  the  strength  of  the  bond  between  the  particles  and  matrix.  Weakly  bonded 
particles  are  pulled  out  in  the  wake  of  the  crack  and  increase  toughness  by  frictional  energy 
dissipation  (McHugh  et  al,  1966;  Hing  and  Groves,  1972).  Particles  firmly  bonded  to  the 
matrix  fracture  in  the  wake  of  the  crack. 

A  number  of  calculations  have  been  done  to  estimate  the  toughening  effect  of  crack 
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bridging.  Many  of  these  axe  concerned  with  frictionaUy  constrained  fibres  in  ceramics 
(Aveston  and  Kelly,  1973;  Hannant  et  al,  1983;  Marshall  and  Cox,  1985;  Budiansky  et 
al,  1986;  Budiansky  and  Amazigo,  1989;  Sigl  and  Evans,  1989).  Other  models  consider 
ceramics  containing  ductile  particles  (Budiansky,  1986,  Budiansky,  Amazigo  and  Evans, 
1988).  Several  studies  are  devoted  to  crack  bridging  in  concrete  (Hillerborg  et  al,  1976; 
Velazco  et  al,  1980;  Wecharatana  and  Shah,  1982;  Ballarini  et  al,  1984).  The  essential 
features  of  these  models  are  similar:  they  consider  a  two-dimensional  (plane  strain)  crack 
with  a  bridged  zone  behind  its  tip.  The  effect  of  this  zone  is  represented  by  a  diet 
of  pressure  p(6),  acting  on  the  crack  flanks.  The  simplest  models  assume  that  p  is  a  linear 
function  of  the  crack  opening  displacements  6  (Rose,  1987).  More  complicated  expressions 
have  been  used  to  model  ductile  particles  or  frictionally  constrained  fibres.  The  two- 
dimensional  calculations  predict  substantial  improvements  in  toughness  due  to  bridging: 
increases  of  10-20  times  the  matrix  toughness  have  been  estimated  for  a  brittle  material 
reinforced  by  ductile  particles,  and  2-4  times  for  a  fibre  reinforced  solid. 

There  are,  however,  some  shortcomings  associated  with  two-dimensional  models  of 
crack  bridging.  The  effect  of  the  crack  front  bowing  between  obstacles  cannot  be  taken 
into  account.  Separate  analyses  suggest  that  crack  trapping  may  considerably  influence 
toughness  (Lange,  1970;  Evans,  1972;  Rose,  1975;  Fares,  1990;  Gao  and  Rice,  1990,  Bower 
and  Ortiz,  1990).  It  is  also  difficult  to  estimate  an  appropriate  form  for  the  p  —  8  relation. 
Particles  properly  bonded  to  the  matrix  are  highly  constrained,  and  so  behave  elastically 
even  when  subjected  to  loads  up  to  6  times  their  yield  stress  (Ashby  et  al,  1989).  The 
particles  therefore  essentially  pin  discrete  points  on  the  faces  of  the  crack.  While  it  is 
sensible  to  represent  the  discrete  pinning  points  as  a  distribution  of  pressure  if  the  length 
of  the  bridging  zone  greatly  exceeds  the  particle  spacing,  this  approximation  breaks  down 
if  there  are  only  a  few  intact  particles  behind  the  crack  front. 

Until  recently,  a  full  three  dimensional  analysis  of  crack  bridging  has  been  beyond  the 
scope  of  conventional  numerical  methods.  However,  progress  in  the  area  of  weight  function 
theory  in  fracture  mechanics  has  provided  an  invaluable  tool  for  solving  problems  of  this 
nature.  Two  different  approaches  have  been  proposed.  The  first  is  a  standard  boundary 
integral  formulation,  where  the  crack  is  represented  by  a  three-dimensional  distribution  of 
dislocations.  The  computing  effort  required  in  this  procedure  may  be  significantly  reduced 
by  using  a  solution  found  by  Rice  (1985)  for  the  interaction  of  a  point  dislocation  with  a 
crack  front  as  a  kernel.  This  technique  was  developed  by  Fares  (1990)  to  solve  the  problem 
of  a  crack  bowing  out  between  obstacles. 

A  second  approach  is  based  on  work  by  Rice  (1985,  1987)  and  Gao  and  Rice  (1986, 
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1987),  who  developed  a  linear  perturbation  scheme  for  calculating  the  variation  in  stress 
intensity  factor  due  to  small  changes  in  crack  geometry.  This  method  was  used  by  Gao  and 
Rice  (1990)  to  find  first  order  solutions  for  the  shape  of  a  crack  front  trapped  by  an  array  of 
tough  obstacles.  Their  analysis  was  restricted  to  the  situation  where  the  particle  toughness 
exceeds  that  of  the  matrix  material  by  a  small  amount,  so  that  the  crack  front  penetrates 
the  obstacles  rather  than  bypassing  them.  More  recently,  Bower  and  Ortiz  (1990)  followed 
a  suggestion  by  Rice  (1987)  to  extend  the  perturbation  method  to  arbitrarily  large  changes 
in  crack  geometry.  They  then  applied  this  technique  to  find  the  shape  of  a  crack  front 
is  it  bypasses  the  array  of  obstacles.  Their  results  were  in  good  agreement  with  those  of 
Far  es  (1990).  In  this  paper,  the  finite  perturbation  method  is  extended  to  take  into  account 
the  effects  of  pinning  particles  in  the  wake  of  the  crack,  and  is  used  in  a  three  dimensional 
analysis  of  crack  bridging. 

The  problem  to  be  solved  is  illustrated  in  Fig  2.  Consider  am  unbounded  elastic  solid 
containing  a  straight,  semi-infinite  crack  on  the  plane  x3  =  0,  oriented  so  that  its  front 
is  parallel  to  the  X2  axis.  The  solid  is  subjected  to  remote  loads,  which  set  up  a  uniform 
distribution  of  mode  I  stress  intensity  K(s)  =  K°°  along  the  crack  front.  Now,  assume 
that  there  is  a  doubly  periodic  distribution  of  particles,  with  spacing  L  in  the  solid  ahead 
of  the  crack  front.  For  simplicity,  the  elastic  constants  of  the  particles  and  matrix  material 
are  taken  to  be  identical,  but  the  fracture  toughness  of  the  particles  A'Par  exceeds  that  of 
the  matrix  A'™at.  The  particles  may  be  cylindrical,  with  the  axis  of  each  cylinder  parallel 
to  the  X3  axis,  this  would  be  representative  of  a  fibre-reinforced  material.  Alternatively, 
the  particles  may  be  spherical,  with  the  centre  of  the  particles  on  the  plane  X3  =  0  (a 
typical  configuration  in  a  ceramic  reinforced  with  metallic  particles).  In  either  case,  the 
r3  =  0  plane  ahead  of  the  crack  front  contains  a  distribution  of  circular  regions,  radius  R, 
where  the  toughness  exceeds  that  of  the  matrix.  Initially,  we  assume  that  the  particles  are 
perfectly  bonded  to  the  matrix  material.  The  analysis  is  subsequently  extended  to  allow 
particles  to  be  pulled  out  in  the  wake  of  the  crack.  We  now  seek  to  calculate  the  behaviour 
of  the  crack  as  it  propagates  through  the  composite. 

Qualitatively,  the  behaviour  of  the  crack  is  as  follows.  When  the  remote  stress  reaches 
a  magnitude  such  that  K°°  =  A™at,  the  crack  begins  to  propagate  through  the  solid. 
Eventually,  parts  of  the  crack  front  run  into  the  tough  particles.  Since  K  \s)  <  Kr , 
these  parts  arrest.  When  the  remote  load  is  increased  by  a  small  amount  the  remainder  of 
the  crack  front  starts  to  bow  out  between  the  pinning  particles,  Fig  1.  The  shape  of  the 
bow  is  determined  by  the  condition  that  K{s)  —  A'™at  over  the  propagating  regions,  while 
K ( s )  <  A'Par  over  regions  in  contact  with  an  obstacle.  There  are  two  possible  ways  in  which 
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the  crack  can  propagate  past  the  row  of  particles.  If  K £ar  is  less  than  a  critical  value,  to  be 
determined,  then  the  particles  are  penetrated  by  the  crack  front  and  eventually  fracture. 
The  sequence  of  trapping,  penetration  and  fracture  repeats  as  the  crack  propagates  through 
the  solid.  Alternatively,  the  crack  may  bow  out  between  the  obstacles  to  such  an  extent 
that  parts  of  its  front  on  the  fax  side  of  the  paticle  begin  to  attract  one  another.  Beyond 
this  point,  the  crack  continues  to  grow  unstably  under  decreasing  remote  load.  Eventually, 
the  crack  joins  up  on  itself,  leaving  behind  a  row  of  circular  bridging  particles.  The  semi- 
infinite  crack  front  then  continues  to  propagate  until  the  next  row  of  obstacles  is  reached, 
whereupon  the  process  is  repeated.  Several  rows  of  bridging  particles  may  form  iii  the  wake 
of  the  crack  in  this  manner.  However,  the  particles  cannot  remain  intact  indefinitely,  and 
may  fail  by  one  of  two  mechanisms.  If  the  particles  are  strongly  bonded  to  the  matrix, 
they  fracture  in  the  wake  of  the  crack  without  appreciable  pull-out.  Alternatively,  if  K £ar 
is  large  compared  to  the  strength  of  the  bond,  the  particles  remain  intact,  and  axe  pulled 
out  in  the  wake  of  the  crack.  If  the  toughness  of  the  interface  is  low,  the  pull-out  is  resisted 
mainly  by  friction  forces  acting  between  the  particle  and  the  matrix,  which  decrease  as  the 
particle  is  pulled  free. 

Analysis  of  this  problem  involves  calculating  stress  intensity  factors  for  a  semi-infinite 
crack  with  an  arbitrarily  shaped  front,  with  its  faces  pinned  over  small  circular  regions. 
It  proves  convenient  to  solve  the  problem  in  stages.  The  solution  begins  with  a  straight 
semi-infinite  crack  subjected  to  remote  loads.  The  bridging  particles  are  then  added,  by 
finding  a  distribution  of  point  loads  on  the  crack  so  as  to  pin  its  faces  over  small  circular 
regions.  This  part  of  the  analysis  is  approximate,  in  that  the  particle  radius  is  assumed 
to  be  much  smaller  than  the  spacing.  Next,  the  crack  front  is  perturbed  from  its  straight 
configuration  using  the  method  described  in  Bower  and  Ortiz  (1990),  until  the  fracture 
criterion  is  satisfied  on  its  front.  Finally,  the  crack  is  marched  around  a  row  of  obstacles  by 
applying  a  further  succession  of  first  order  perturbations.  In  this  way,  a  complete  picture 
of  the  behaviour  of  the  crack  as  it  bypasses  each  row  of  obstacles  can  be  built  up. 


2.  Theory 

2.1  Straight  crack  perfectly  pinned  by  arrays  of  particles 


Consider  the  straight,  semi-infinite  crack  shown  in  Fig.  2  subjected  to  remote  loading, 
so  that  a  uniform  mode  I  stress  intensity  factor  K(s)  =  K°°  is  induced  on  the  crack  front. 
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Initially,  the  crack  is  free  of  pinning  particles,  and  the  displacement  at  x  is  A U(x),  with 

(2i 

A  _  (1  _  ^ )  (0  r. 


where  v  is  Poisson’s  ratio  and  p  the  shear  modulus.  Now,  a  single  bridging  particle  may  be 
introduced  in  the  wake  of  the  crack  by  finding  an  appropriate  distribution  of  point  loads 
over  the  pinned  region.  First,  note  that  the  displacement  at  x  due  to  a  pair  of  unit  wedging 
forces  acting  at  £  is  A G(£,x)  (Ulfyand,  1965),  where 


x)  =  ^ arctan j2 y  ~y  f  1  (2-3) 

P  =  [(^i  —  £i)2  +  (*2  —  £2 )2]  ^  • 

Distributing  wedging  forces  over  a  circular  domain  ft,  radius  R,  centered  at  x^  on  the 
crack  faces  then  gives  rise  to  the  following  displacement  at  x 


v(x)  =  AU(x)  +  A  f  G(S,x')p(Z)dt. 
J  n 


If  the  crack  faces  are  completely  pinned  by  the  bridging  particle,  then  v(x)  —  0  for  x  €  ft, 
and  (2.4)  gives  an  integral  equation  for  the  unknown  distribution  of  wedging  forces  p(£)  on 
the  particle.  Note  that  it  is  not  necessary  to  enforce  complete  closure  over  ft:  any  arbitrary 
relation  between  v  and  p  could  be  prescribed.  In  this  way,  the  analysis  may  be  extended 
to  model  debonding  between  the  particle  and  the  matrix,  or  plastic  deformation  in  the 
particle.  This  issue  is  addressed  in  section  2.4,  but  for  simplicity  we  begin  by  assuming 
that  the  particle  is  perfetly  bonded  and  u(x)  =  0. 

Equation  (2.4)  cannot  be  solved  analytically  in  its  present  form.  However,  if  the  par¬ 
ticle  radius  is  much  smaller  than  the  distance  to  the  crack  front,  G(£,x)  may  be  simplified 
as  follows.  Begin  by  expressing  G  in  terms  of  polar  co-ordinates  centered  at  x(n),  Fig.  2 


G(r',6',  r,  9)  =  — ^  arctan  <  2\ 


<  (xY  —  r  cos  Q)(x\n ^  ~  r' cos  &') 


p  —  [r2  4-  r'2  —  2 rr'  cos (9  —  6' )]  ^  . 

Expanding  the  arctan  in  (2.5)  to  first  order  in  p/x\n\  and  letting  r/x[n\  r'  jx^  — *  0  then 


leads  to  the  following  approximation  to  G 


G(i,x) 


1  (  K 

"T-  «  + 


\  “  2x 
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Similarly,  we  may  write  £7(x)  ~  U(x<'n 1)  to  the  same  order  of  accuracy.  The  first  term  on 
the  right  hand  side  of  (2.6)  represents  the  displacement  field  in  a  half-space  due  to  point 
loading.  The  second  term  corrects  the  solution  for  the  presence  of  undamaged  material 
ahead  of  the  crack  front. 

For  R  <<  the  integral  equation  (2.4)  therefore  reduces  to 


0  =  P(x(n))  + 


1  /  (- 


+ 


J n  \2 p  2x\ 


(") 


p(p,  6)  dpd.6. 


(2.7) 


A  similar  equation  arises  in  the  analysis  of  a  rigid  circular  punch  indenting  a  half-space, 
with  p(£)  as  the  unknown  pressure  between  the  contacting  surfaces,  and  At/(x^)  their 
normal  approach  (Ulfyand,  1965).  The  solution  has  the  form 


p(p )  = 


47tP2 


R  +  P  R  ~  P 
R  —  p  V  R  +  p 


(2.8) 


P  = 


(xl,”l-6)2+(4",-6) 


(*») 


Sen. 


where  P  is  the  resultant  of  the  pressure  distribution.  Substituting  for  p(p)  in  (2.7)  and 
integrating,  we  obtain  an  expression  which  may  be  solved  for  P 


0  =  P(x(n))  + 


P  P  \ 
«  +  2^x[n) )  ' 


(2.9) 


Here,  P/4P  represents  the  displacement  under  a  circular  punch  indenting  a  half-space  and 
is  independent  of  the  position  of  the  crack.  The  second  term  corrects  the  solution  for  the 
constraining  effect  of  the  unbroken  surface  ahead  of  the  crack. 

Now,  noting  that  the  pressure  distribution  (2.8)  is  square  root  singular  at  p  =  R,  the 
stress  intensity  on  the  crack  front  mound  the  pinning  particle  may  be  calculated  in  the 
usual  way,  with  the  result 

Arpar  = - (2.10) 

2R\ZttR 

The  negative  sign  in  equation  (2.10)  is  a  consequence  of  the  sign  convention  for  G.  which 
was  chosen  so  that  the  positive  forces  acting  on  the  crack  faces  tend  to  open  the  crack.  The 
pressure  on  the  pinning  particles  is  compressive,  so  P  <  0.  Equation  (2.10)  corresponds  to 
the  well  known  solution  for  the  stress  intensity  on  an  external  circular  crack  subjected  to 
uniform  remote  tension. 


Finally,  the  stress  intensity  factor  on  the  straight  crack  front  may  be  corrected  for 
the  presence  of  the  pinning  particle.  To  first  order  in  R/x[n\  the  pinning  particle  may  be 


9 


approximated  by  a  pair  of  wedging  forces,  magnitude  P,  acting  on  the  crack  faces.  They 
induce  a  stress  intensity  factor  K(s)  —  PH(s,x(n 1),  where 


H(s,x(n))  = 


^  (52+(x^)2) 


(2.11) 


The  total  stress  intensity  factor  is  therefore  given  by 


K(s)  =  K°°  +  PH(s,x(n >). 


(2.12) 


It  is  straightforward  to  extend  the  preceding  analysis  to  allow  for  an  arbitrary  number 
of  pinning  particles:  the  configuration  shown  on  Fig.  2  is  of  particular  interest.  The  crack 
is  now  pinned  by  N  parallel  rows  of  circular  particles.  The  spacing  between  particles  in  the 
i2  direction  is  L,  but  otherwise  their  positions  are  arbitrary,  and  may  be  specified  in  terms 
of  the  co-ordinates  x^"'  of  a  chosen  particle  on  the  nth  row. 

We  may  write  down  a  closure  condition  analogous  to  equation  (2.9)  for  each  row  of 
particles.  It  is  only  necessary  to  add  the  effects  of  interactions  between  the  particles.  To 
this  end,  we  make  use  of  the  additional  assumption  that  the  spacing  between  particles 
greatly  exceeds  their  radius,  so  that  the  displacement  at  x  due  to  a  particle  at  £  may  be 
found  by  replacing  the  particle  with  a  pair  of  wedging  forces.  This  immediately  leads  to 


(1  " 

—  + - — +  Y,  G(x(n>+mXj,x(n)) 


1  TTI  —  —  OO 

mqtO 


+  Eptk)  E  G(xw  +  ml,x(n)), 


(2.13) 


k  =  l  T/i—  ■*“  O./ 

k^n 


where  P^  is  now  the  pinning  force  on  a  particle  in  the  nth  row,  and  j  is  a  unit  vector 
in  the  r2  direction.  It  will  prove  convenient  later  to  simplify  this  result  by  defining  a  new 
function  G  given  by 


G(x("\x(fc))  = 


Y  G’(x(i)  +  mlj,x(n)), 


n  ^  k. 


(2.14) 


Clearly,  the  function  G(x^fc\x^nl)  represents  the  displacement  at  a  particle  in  the  the  nth 
row  due  to  the  pinning  forces  acting  on  the  fcth  row.  Equation  (2.131  may  now  be  written 

N 

0  =  17(x(n))  +  ^P(*)G(x(n),x(t)).  (2.15) 

*=i 
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giving  a  set  of  N  linear  equations  which  may  be  solved  for  the  unknown  P(n>.  The  total 
number  of  rows  of  particles  N  may  be  found  from  the  condition  that  the  stress  intensity 
factor  on  the  last  row  must  not  exceed  the  particle  toughness.  From  equation  (2.10),  this 


\PiS)\  <  2RV^RKTt. 


(2. 10) 


Similarly,  the  expression  for  the  corrected  stress  intensity  factor  on  the  straight  crack 
front  (equation  (2.12))  becomes 

N  +no 

K{s)  =  K°°  +  Y  P{n)  Y  H(s,x{n)  +mlj),  (2.17) 


n=  1  m=  — oc 


which  may  be  simplified  by  putting 


HLs.x(n))=  Y  H(s,x{n)  +  mlj)  = 


-2ri  (sinh )/(2~~.  i 
~3  (£/;r)2  (sixth2  4-  sin*  ~rj) 


/here  *•  =  —i\n)/L  and  rj  -  (s  -  r!/*1  )j  L  to  give 

N 

K(s)  -  A' 30  +  ^P(n)H(s.x(’‘) 


It  is  now  possible  to  estimate  the  magnitude  of  the  remote  load  (paramaterised  by 
A’°°)  required  to  advance  the  crack  through  the  material.  Strictly  speaking,  it  is  necessary 
to  enforce  the  condition  I\(s)  —  A'™at  at  all  points  on  the  crack  front  at  fracture.  However, 
once  the  bridging  particles  have  been  added  to  the  wake  of  the  crack,  the  c  tress  intensity 
factor  on  its  front  is  no  longer  uniform.  Consequently,  it  is  not  possible  to  satisfy  the 
fracture  criterion  exactly  without  changing  the  shape  of  the  crack  front.  This  issue  is 
addressed  in  the  next  section.  Nevertheless,  it  is  possible  to  obtain  an  accurate  estimate 
of  the  increase  in  toughness  due  to  bridging  by  setting  the  average  stress  intensity  factor 
along  the  crack  front  equal  to  A”™at.  It  is  straightforward  to  show  that  the  mean  value  of 
K(s)  is: 

1  P(n)  /  _ 9 

<*<*)>- z/  +  ,/^T  <22°> 

n  =  1 


(2.20) 


whence  the  toughening  ratio  given  by 


mat 

c 


1  + V^E— p~t= 

V  n  ^  T  T'r /  (n) 
n=l  LI\  00  \ j  X\ 
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By  way  of  illustration,  it  may  be  shown  that  for  the  special  case  of  a  single  row  of 
particles,  distance  b  behind  the  crack  front 


where 


pil)  _  /26  f  L _ L_  2F(b/L)  | 

LK°°  V  7T  [  4R  2n 2b  7r  J 

A' 00  _  irL/(8R)  -  I/(4tt6)  +  F(b/L) 

A'mar  ~  nL/{8R)  -  A/(4r6)  +  F(b/L)  -  r 

A'par  _  L  nr  1 

~  2R\  2R  7rL/(8A)  -  L/(4tt{>)  +  F(b/I)  -  1  ’ 


no  -  E 


1 


m=  1 


1 

—  arctan 
m 


(2.22) 

(2.23) 

(2.24) 

(2.25) 


2.2)  Extension  to  wavy  cracks 

It  is  evident  from  the  preceding  discussion  that  the  crack  front  will  have  a  wavy  profile  at 
fracture.  In  addition,  parts  of  the  crack  may  bow  out  between  trapping  particles.  In  order 
to  take  into  account  these  effects,  it  is  necessary  to  extend  the  analysis  to  cracks  with  an 
arbitrarily  shaped  front:  this  issue  is  addressed  next. 

The  following  approach  will  be  used.  The  analysis  described  in  the  previous  section 
makes  use  of  three  functions:  A(x).  giving  the  displacement  at  x  under  a  uniform  remote 
load;  G(x,£) ,  which  gives  the  displacement  at  x  due  to  point  loads  applied  at  and 
H(s.x).  which  gives  the  stress  intensity  factor  at  a  point  s  on  the  crack  front  due  to 
wedging  forces  at  x.  Clearly,  if  U  G  and  H  could  be  found  for  an  arbitrarily  shaped 
crack,  the  procedure  outlined  in  section  2.1  could  be  applied  in  the  same  way  to  add  rows 
of  pinning  particles.  In  practice,  the  functions  are  known  exactly  only  for  circular  and 
semi-infinite  cracks.  However,  it  is  possible  to  calculate  the  change  in  the  functions  caused 
by  a  small  perturbation  to  the  geometry  of  the  crack,  using  the  procedure  developed  by 
Rice  (1985).  By  applying  a  succession  of  such  perturbations,  G ,  H  and  U  may  then  be 
progressively  updated  for  a  crack  with  an  arbitrarily  shaped  front.  The  equations  involved, 
cannot  be  solved  analytically,  so  a  numerical  procedure  is  used. 

We  begin  by  summarising  Rice's  equations  for  the  effects  of  a  first-order  perturbation 
in  geometry.  For  details  of  the  derivation  of  these  results,  the  reader  is  referred  to  the  work 
of  Rice  (1985,1987),  Gao  and  Rice  (1986,1987)  and  Bower  and  Ortiz  (1990).  Consider  a 
planar  semi-infinite  crack  with  an  arbitrarily  shaped  front,  and  assume  that  the  functions 
G(x,  £),  H(s,£)  and  U(£)  are  known  for  this  configuration.  In  addition,  assume  that  the 
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ciack  is  subjected  to  remote  loading,  so  that  there  is  an  initial  mode  I  stress  intensity  factor 
K(s)  on  the  crack  front.  The  magnitude  of  the  load  will  be  parameterised  by  a  remotely 
applied  stress  intensity  factor  A'00,  defined  by  considering  the  energy  released  during  an 
infinitessimal  crack  advance  in  the  xj  direction.  Given  K(s),  K°°  is  found  to  be 


K°°  =  y  —  J0  K2^dl2'  (2'2G) 

where  L  is  one  wavelength  of  the  crack  front. 

It  is  convenient  at  this  point  to  introduce  another  weight  function  for  the  crack,  defined 
by 

_  fn*  H( s.z(t)  —  en{i)) 

D(s,t)  =  limt/  — — - - ’-L,  (2.27  i 

^-OV  2  yft 

where  z (t)  is  the  position  vector  of  a  point  t  on  the  crack  front  C,  and  n (t)  is  the  normal  to 
C.  Evidently.  D  represents  the  stress  intensity  factors  at  a  point  s  on  the  crack  front  as  a 
pair  of  square  -root  singular  wedging  forces  approach  the  crack  front  at  t.  The  significance 
of  this  weight  function  is  discussed  in  Bower  and  Ortiz  (1990).  For  instance,  using  (2.1i ). 
it  is  easy  to  show  that  for  a  straight  crack 


D(s,t) 


{t  -  s) 


Now.  imagine  allowing  the  crack  to  advance  by  a  small  amount  <5a(s)n(s).  The  follow¬ 
ing  equations  may  then  be  shown  to  hold,  accurate  to  first  order  in  Sa(s): 

SK00  1  f 

6U(x)  =  +  -  /  K(s)H(s,x)6a(s)  ds.  (2.29) 


£K° c  1  f 

6U(x)  =  -j-r^-U(x)  +  -  /  K(s)H(s,x)6a(s)ds. 

J  c 

6G(x,£)=-  I  H(s,£)H(s,x)6a(s)  ds, 

2  Jc 

6K°°  1  f 

6K(s)  —  —  -  K0(s)  +  —  /  K{t)D{t,s)[6a(t)  —  n{s)  •  n{t)6a{s)]dt, 

A  00  2.n  Jc 

6H(s,  x)  =  —  J  H(\,x)D{s,\)^6a{\)~ 


[x  -  z(A)]  •  [x  -  z(s)] 
[x  -  z(s)]  •  [x  -  z(s)] 


n(A)  ■  n(s)(5a(s)  >  d A 


Sa(s) 

2tt 


f  E7 /  \  rn/  ,  J(x -z(A))  x  (x  -  z(s))]  ■  kr  ,  „  ,  JX 

Jc  H(\,  x)D(s.  A)  z(  -  -  [n(A)  x  „(*)]  -  kd\ 


2t t  Jc  '  '  [x  -  z(s)J  •  [x  -  z(s)] 

3 Sa(s)  [x  -  z(.s)]  •  n(s) 

6D(s.t)  =  ~  J  D(A.t)D(s,A)|M*!- 


(2.30) 


(2.31) 


(2.32) 
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[*(*)  -  *(0J  •  [2(A)  -  *(*)] 


6a{s)  J^D(X,t)D(s,X) 

/  D(X,t)D(s,X) 
Jc 


2  7T 
6a(t) 
2tt 


-  2 D(s,t) 


[z(0  -  z(s)]  •  [z(t)  -  z{s)} 

_  (g(0  ~  z(3)l  •  iz(A)  ~  z(5)l 

[z{t)  -  z(s)]  •  [z(<)  -  z(«)] 
[z(A)  -z(f)]  x  [z (s)  -  z(t)]  •  k 
[z(t)  -  z(j)j  •  [z (t)  -  z(s)] 
(z(A)  -z(a)]  x  [z(t)-z(j)]'k 
[z(s)  -  z (t)}  •  (z(s)  -  z(t)] 
[£a(f)n(f)  —  6a(s)n(s)]  •  [z(t)  —  z(s)| 


n(A)  •  n (s)<5a(s) 

n(A)  •  n(.s)<5a(f)| 


dX 


n(A)  x  n(s)  •  kdX 
n(A)  x  n(f)  •  kdX 


[z (0  -  z (s)}  ■  [z (t)  -  Z(5)] 


(2.33) 


Here,  fc  denotes  integration  around  the  crack  front;  (x)  and  (•)  are  vector  and  scalar 
products.  In  equations  (2.29)  and  (2.31),  the  terms  involving  6K°° / K°°  represent  changes 
due  to  an  increment  ;n  the  remote  load;  the  remaining  terms  are  changes  due  to  the 
perturbation  in  crack  geometry.  Equations  (2.32)  and  (2.33),  which  are  modified  forms 
of  Rice’s  results,  are  derived  in  Appendix  A.  The  remaining  equations  may  be  found  in 
Rice  (1985,  1987). 

These  results  are  the  basis  of  an  incremental  method  for  calculating  stress  intensity 
factors  and  weight  functions  for  any  arbitrarily  shaped  crack.  For  example,  stress  intensity 
factors  may  be  calculated  by  repeatedly  applying  equations  (2.31)  and  (2.33).  The  analysis 
begins  with  a  convenient  reference  geometry  for  which  stress  intensities  and  influence  func¬ 
tions  are  known  in  closed  form,  such  as  a  circular  or  semi-infinite  crack.  Equation  (2.31)  is 
then  used  to  calculate  the  change  in  stress  intensity  caused  by  a  small  change  in  the  shape 
of  the  crack.  The  problem  may  be  posed  in  one  of  two  ways:  we  may  either  find  changes 
due  to  a  prescribed  crack  advance;  or  alternatively,  we  may  calculate  the  displacement  of 
the  crack  front  consistent  with  a  given  crack  growth  criterion,  such  as  K(s)  =  Kc(z(s)). 
In  the  latter  case,  (2.31)  is  a  singular  integral  equation  to  be  solved  for  6a(s).  Once  the 
new  shape  of  the  crack  has  been  found,  equation  (2.33)  is  used  to  calculate  a  new  weight 
function  D  for  the  perturbed  geometry.  Similarly,  the  crack  face  displacements  U ,  and  the 
weight  functions  G  and  H  may  be  updated  using  equations  (2.29),  (2.30)  and  (2.32).  This 
process  may  be  repeated  ad  infinitum:  by  successively  updating  the  crack  geometry  and 
influence  functions,  we  may  calculate  the  effect  of  an  arbitrarily  large  displacement  of  the 
crack  front.  In  actual  calculations,  the  governing  integral  equations  are  discretised  using 
standard  numerical  techniques.  A  method  of  solution  has  been  described  by  Bower  and 
Ortiz  (1990),  and  is  summarised  briefly  in  Section  2.5.  Since  the  appropriate  equations  are 
accurate  to  first  order  in  <5a(s),  the  numerical  solution  may  be  expected  to  converge  as  the 
amplitude  of  each  successive  perturbation  and  the  mesh  size  are  progressively  reduced. 
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For  the  particular  case  of  a  periodic  perturbation  of  a  semi-infinite  crack,  equations 
(2.29)  to  (2.31)  may  be  simplified  by  reducing  the  integrals  to  one  wavelength  of  the  crack 
front.  In  addition,  it  is  possible  to  derive  incremented  expressions  for  the  summed  weight 
functions  G  and  H  defined  in  equations  (2.14)  and  (2.18).  The  derivations  are  somewhat 
laborious,  and  are  given  in  appendix  B.  Here,  we  merely  summarise  the  results.  Consider 
allowing  the  front  of  a  semi-infinite  crack  to  advance  by  a  small  amount  Sa(s),  which  is 
periodic  in  the  xi  direction  so  that 


6a(s  +  me )  =  6a(s),  (2.34) 


where  c  is  the  arc  length  along  one  wavelength  (c  =  L  for  a  straight  crack),  and  m  is  any 
integer.  It  then  follows  that 

SU(x)  =  6J£-U(x)+  J  f/2  K0(s)H(s,x)6a(s)ds,  (2.35) 

K  Z  J-c/2 

1  rcl 2 

6G{x{n\x(k))  =  -  H(s,x(n))H(s,  x^k^)6a(s)  ds,  (2.36 ) 

2  J-c/2 

sk°°  i  rc /2 

SK(s)  =  -bsr- *(s)  +  7T  /  K(X)D(\,  i)NA)  -  6a(s) n(A)  •  n (s)}  dX,  (2.37) 

^  J-c/2 

1  fc/ 2 

6H(s,x)  =  -^  /  H(A,x)D(A,  s)[<5a(A)  —  <5a(s)n(A)  •  n(s)]  dX 

2?r  J —c/2 

+6a(s)'&(s,x),  (2.38) 


where: 


+  oo 


D ($,t)  =  E  D(s  4-  me,  t) 


i  pc/2  +  +°° 

^'X)  =  s/.„E  E  +  (m  —  n)Lj)L>(A  +  nc,  s) 

J  ~c/2  n  =  —  oc  m=  —  txj 

z(s)  -  z(A)  +  -nLj]  •  [x  -  z(s)  +  mLj] 


[x  -  z(s)  +  mLj]  ■  [x  -  z(s)  +  mLj] 

[x  -  z(A)  4-  (m  —  n)Lj]  x  [x  —  z(s)  4-  mLj]  •  k 


[x  -  z(s)  +  mLj]  •  [x  -  z (s)  4-  mLj] 


n(A)  •  n (s)} 

[n(A)  x  n(s)]  •  k  J 


dX 


J  o  ,  r;,  [x-z(5)4-mLj]-n(a) 

?  2_^  s,  x  m  +  mLj]  •  [x  -  z(s)  4-  mLj] 


(2.39) 


(2.40) 


The  updated  influence  functions  G  and  H  may  be  used  to  add  pinning  particles  to 
the  wake  of  the  crack  as  outlined  in  section  2.1  For  a  semi-infinite  crack  pinned  by  N  rows 
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of  parallel  circular  particles,  the  pinning  force  on  each  row  of  obstacles  is  calculated  using 
the  closure  condition  (2.15).  Here,  we  have  assumed  without  proof  that  for  R/L  <<  1,  the 
distribution  of  pressure  on  each  pinned  region  is  unchanged  by  modifying  the  crack  front, 
and  may  still  be  characterised  by  its  resultant  P.  This  is  justified  in  Appendix  B.  The 
stress  intensity  factors  on  the  crack  front  are  then  found  from  (2.19) 

In  general,  once  pinning  particles  have  been  added  to  the  wake  of  the  crack,  it  is 
necessary  to  modify  the  crack  geometry  in  order  to  satisfy  the  fracture  criterion  at  all 
points  on  the  crack  front.  For  example,  it  was  shown  in  section  2.1  that  when  pinning 
particles  are  added  to  a  straight  crack,  the  condition  K{s)  =  A'“at  can  only  be  met  by 
perturbing  the  crack  front.  This  issue  is  addressed  next. 

The  problem  may  be  stated  as  follows.  We  wish  to  find  a  crack  advance  8a(s)  that 
changes  the  stress  intensities  on  the  crack  front  by  a  known  amount  6K(s),  such  that 
K(s)  -f  6K(s)  =  Kc(z(s)).  At  the  same  time,  we  wish  to  satisfy  the  closure  condition  on 
each  particle  for  the  perturbed  geometry.  To  this  end,  we  write  down  expressions  similar 
to  (2.15)  and  (2.19)  for  the  perturbed  crack: 

N 

0  =  U(x(n])  +  8U(x(n))  +  £  ( P(k)  +  6Pik)^j  (G(x<*\x<n>)  +  <5G(x(i),x(n)))  ,(2.41) 

N 

K(s)  +  6K{$)  =  K0(s)  +  6K0(s)  +  ( P <*>  +  <5P(*>)  (H(5,x<n>)  +  <5H(s,x(n>));?.42) 

where  6P(*)  is  the  change  in  pinning  force  on  the  fcth  row  of  particles.  Neglecting  second 
order  terms,  and  using  (2.15)  and  (2.19),  these  results  may  be  simplified  to 

N 

0  =  6U(x(n) )  +  ^2  (6p(k) G(x(n),x(*>)  +  P{k)6G(x^\x{k)))  ,  (2.43) 

k=l 

N 

6K(s)  =  6K0(s)  +  ]T  {6p(k) H(^x(n))  +  PC*^H(s,x(n)))  .  (2.44) 

*=i 

Finally,  substituting  (2.35)  to  (2.38)  into  (2.43)  and  (2.44),  and  eliminating  Ko(s)  and  U(x) 
using  (2.15)  and  (2.19),  we  find  that 

0  =\J//  K(\)H(\,x(n))6a(\)d\  +  J2  ~  G(x(n\x(t)),(2.45) 

1  /"C/2  £  l(OC 

6K(s)  =  —  /  K(X)D(X,s)[6a(X)  -  n(A)  •  n(s)6a(s)]  dX  + 

a7T  J  —  c/2  A 

N  N  /  c  jsoo  \ 

+  6a(s)  ^  P{k)V(s,x{k})  +  E  6p(k)  ~  ~F^p{k)  )  H(5,x^>).  (2.46) 

*=i  )t=i  '  \°°  J 
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Here,  K(s)  represents  the  stress  intensity  factor  on  the  pinned  crack:  all  variables  relating 
to  the  unpinned  crack  have  been  eliminated.  Given  an  increment  in  stress  intensity  6K(s), 
(2.45)  and  (2.46)  form  a  set  of  N  4- 1  coupled  singular  integral  equations  to  be  solved  for 
6a(s)  and  <5PO).  Once  the  new  shape  of  the  crack  has  been  found,  equations  (2.32),  (2.33), 
(2.36)  and  (2.38)  may  be  used  to  calculate  updated  influence  functions  for  the  perturbed 
geometry. 


2.3)  Application  to  crack  growth  through  toughened  solid 


The  procedure  outlined  in  the  previous  section  may  now  be  applied  to  solve  the  specific 
problem  of  a  semi-infinite  crack  propagating  through  a  material  containing  tough  second 
phase  particles.  Consider  the  configuration  shown  in  Fig.  2.  A  particulate  toughened  solid 
contains  a  straight,  semi-infinite  crack:  initially,  the  crack  faces  are  free  of  pinning  particles, 
and  we  will  assume  that  the  crack  front  is  just  short  of  one  row  of  obstacles.  Now,  imagine 
applying  a  steadily  increasing  remote  load  to  the  body.  When  the  load  reaches  a  magnitude 
K°°  =  Ar™at,  the  crack  starts  to  propagate  around  the  particle.  Parts  of  the  crack  front 
which  touch  the  particle  are  assumed  to  arrest,  so  that  6a(s)  =  0  at  these  points.  The 
crack  advance  over  the  remainder  of  the  crack  front  is  determined  by  the  condition  that 
K(s)  =  A'™at:  this  may  be  satisfied  by  solving  equation  (2.37)  with  6K(s )  =  0.  If  there 
are  no  pinning  particles  on  the  crack  faces,  we  obtain  the  following  conditions  on  the  crack 
advance 


<5a(s)  =  0 


0  = 


6 1< 30 

K°° 


K(s) 


Crack  trapped, 


1  f c/ 2 

+  —  /  K(\)D(\,s)[6a{\)  -  n(\)  ■  n(s)6a(s)]d\ 

27r  J —c/2 


(2.47) 

Crack  propagating.  (2.48) 


In  principle,  given  an  increment  in  remote  load  6K°°,  (2.47)  and  (2.48)  may  be  solved 
for  ^a(s).  However,  this  may  be  an  ill-conditioned  problem,  since  a  very  small  change 
in  load  may  cause  a  large  crack  advance,  particularly  when  the  peak  load  is  approached. 
In  practice,  it  is  preferable  to  leave  the  increment  in  load  as  an  unknown,  and  enforce 
an  additional  constraint  on  the  crack  advance.  Here,  we  have  chosen  to  constrain  the 
maximum  crack  front  displacement.  The  numerical  procedure  used  to  solve  the  singular 
integral  equation  (2.48)  is  described  in  Bower  and  Ortiz  (1990),  and  is  summarised  briefly 
in  the  next  section. 


Once  6a(s)  has  been  determined  over  the  propagating  regions  of  the  crack,  8K(s) 
on  the  arrested  parts  may  be  found  using  (2.37),  and  the  weight  function  D  is  updated 
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using  (2.33).  In  addition,  the  crack  face  displacements  Z/(x),  and  the  influence  functions 
G(x,  £)  and  H(x,s )  are  updated  at  chosen  points  to  allow  pinning  particles  to  be  added 
to  the  crack  at  a  later  stage.  This  procedure  is  repeated  to  march  the  crack  front  around 
the  first  row  of  obstacles.  The  behaviour  of  the  crack  as  it  bypasses  a  row  of  particles  is 
shown  in  Fig.  4.  As  the  crack  front  bows  between  the  inclusions,  it  eventually  reaches  a 
configuration  where  parts  of  the  crack  begin  to  attract  one  another.  At  this  point,  the 
applied  load  reaches  a  maximum,  and  the  crack  subsequently  propagates  under  decreasing 
remote  load.  Eventually,  the  crack  joins  up  on  itself  on  the  fax  side  of  the  obstacle,  leaving 
behind  a  bridging  particle. 

It  is  not  possible  to  analyse  the  coalescence  of  the  crack  using  the  finite  perturbation 
method.  The  perturbation  procedure  only  gives  accurate  results  if  the  change  in  weight 
function  SD  is  small  compared  to  D  during  each  step.  Equation  (2.33)  shows  that  SD 
varies  as  (fz(s)  —  z(t)]  •  [z(s)  —  z(t)])-1,  where  z(s)  and  z (t)  are  the  position  vectors  of 
points  t  and  s  on  the  crack  front.  If  the  crack  is  about  to  join  up  on  itself,  z(s)  — ►  z it)  at 
the  points  of  coalescence.  Consequently,  in  order  for  SD  to  be  small,  <5a(A)  must  become 
vanishingly  small  at  these  points,  and  an  infinite  number  of  time  steps  are  required  for  the 
crack  to  join  up.  In  addition,  D(s,t)  becomes  badly  behaved  at  the  point  of  coalescence, 
so  an  extremely  fine  crack  front  mesh  is  required  to  get  accurate  numerical  results,  making 
the  perturbation  method  prohibitively  expensive. 

These  difficulties  may  be  avoided  by  using  the  method  outlined  in  the  section  2.1  to 
add  pinning  particles  to  the  wake  of  the  crack.  We  therefore  ignore  the  details  of  the 
coalescence,  and  consider  the  crack  at  a  slightly  later  stage.  For  example,  we  may  resume 
the  analysis  at  the  point  when  the  crack  front  just  short  of  the  second  row  of  obstacles,  with 
a  single  row  of  bridging  particles  in  its  wake.  A  crack  configuration  must  now  be  found 
that  satisfies  the  appropriate  displacement  boundary  conditions  on  each  pinning  particle, 
and  the  fracture  criterion  K{s )  =  K™*1  on  the  crack  front.  This  may  be  reached  in  three 
stages: 

(i)  Starting  with  a  straight  crack  in  the  required  position  under  unit  remote  load,  equation 
(2.15)  is  solved  to  find  the  magnitude  of  the  pinning  force  on  each  particle; 

(ii)  The  magnitude  of  the  remote  load  K°°  is  adjusted  so  that  the  fracture  criterion  K(s)  = 
K ™al  at  some  arbitrary  point  (say  s  =  ±c/ 2). 

(iii)  Finally,  the  crack  front  is  allowed  to  advance  until  the  fracture  criterion  is  satisfied  at  all 
points  its  front.  To  do  this,  we  set  SK(s)  =  K(s)--K^'a,t  in  (2.46),  and  solve  the  coupled 
equations  (2.45)  and  (2.46)  for  6a(s ),  subject  to  6K°°  ~  0.  If  the  necessary  SK(s)  is 
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large,  several  small  perturbations  may  be  applied  to  reach  the  correct  configuration. 

This  procedure  may  be  applied  to  add  pinning  particles  to  the  crack  at  any  stage  as 
it  bows  out  between  a  row  of  trapping  obstacles.  A  complete  picture  may  then  be  built  up 
as  the  crack  propagates  around  several  successive  rows  of  obstacles.  This  approach  signifi¬ 
cantly  reduces  the  computing  effort  involved  in  the  analysis,  since  any  arbitrary  number  of 
rows  of  obstacles  may  be  added  without  having  to  solve  for  all  intermediate  configurations 
of  the  crack.  The  actual  number  of  rows  of  obstacles  left  intact  in  the  wake  of  the  crack  is 
determined  by  the  particle  toughness.  Particles  are  added  until  the  stress  intensity  factor 
on  the  particles  furthest  from  the  crack  front  (calculated  using  equation  (2.10))  exceeds 
Kf“. 

2.4)  Frictional  crack  bridging 

In  the  preceding  sections,  it  was  assumed  that  the  particles  remain  perfectly  bonded  to 
the  matrix  material,  and  fail  by  fracture  in  the  wake  of  the  crack.  If  the  interface  between 
particles  and  matrix  is  weak,  it  may  fail  before  the  particles  fracture.  The  particles  then 
remain  intact  and  are  pulled  out  in  the  wake  of  the  crack.  The  pull-out  is  resisted  by 
friction  forces  acting  between  the  particles  and  matrix,  which  decrease  as  the  particle  is 
pulled  free.  The  energy  dissipated  as  long  grains  are  pulled  out  in  the  wake  of  a  crack  is 
believed  to  be  an  important  toughening  mechanism  in  monolithic  ceramics  (Vekinis  et  al, 
1990).  In  this  section,  the  analysis  of  a  pinned  crack  is  extended  to  account  approximately 
for  frictional  pull-out. 

The  separation  of  an  individual  grain  from  the  surrounding  material  is  idealised  using 
the  simple  model  illustrated  in  Fig.  3.  The  grain  is  taken  to  be  cylindrical,  radius  R  and 
height  <5crit.  We  assume  that  the  particle  is  pinched  into  the  matrix  by  thermal  stresses,  so 
that  a  local  compressive  stress  crr  acts  perpendicular  to  the  interface.  The  force  required 
to  pull  the  grain  free  by  a  distance  6  is  given  by 

P  =  p0(  1  _  J-).  (2.49) 

Vcrit 

Here,  Po  is  the  maximum  pinning  force  on  the  grain,  given  by 

P0  =  TrforR6CIit,  (2.50) 

where  /  is  the  coefficient  of  friction.  We  note  in  passing  that  it  is  possible  to  enforce  any 
non-linear  relation  between  crack  opening  and  pinning  force.  For  example,  the  analysis 
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could  be  extended  to  model  plastic  deformation  in  the  particles.  However,  this  is  not 
pursued  here. 

Suppose  that  the  matrix  contains  a  regular  distribution  of  similar  grains,  spacing  L. 
For  simplicity,  assume  that  the  grains  are  all  aligned,  with  axes  parallel  to  the  x$  direction. 
Now,  consider  the  behaviour  of  a  semi-infinite  crack  in  the  X3  =  0  plane  as  it  propagates 
through  the  solid.  The  particles  trap  the  front  of  the  crack  as  before,  which  bows  out 
between  them.  We  assume  that  the  particles  are  sufficiently  tough  never  to  be  penetrated 
by  the  crack,  so  that  they  are  left  behind  in  its  wake.  The  particles  remain  perfectly  pinned 
until  the  pinning  force  exceeds  Po,  at  which  point  they  start  to  be  separated  from  the 
matrix.  During  pull-out,  the  pinning  force  P  is  related  to  the  crack  opening  displacement 
<5  by  equation  (2.49).  When  6  >  <5Crit,  the  particle  is  completely  free  from  the  matrix  and 
P  =  0. 

It  is  straightforward  to  incorporate  this  force-crank  opening  displacement  relation  into 
the  model  of  a  bridged  crack.  By  relaxing  the  constraint  that  u(x)  =  0  in  equation  (2.4), 
the  separation  of  the  crack  faces  is  related  to  the  pinning  force  by 

N 

<5(x(n))  =  2At/(x(n))  +  2A^P(fc)G(x(n),x(fc)).  (2.51) 

jfc=i 

Here  6(x^)  is  the  crack  opening  displacement  at  the  nth  row  of  particles,  A U(x^n^)  is  the 
displacement  at  the  same  point  for  the  unpinned  crack,  and  AG(x^n\  x^)  represents  the 
displacent  at  x^  due  to  the  pinning  forces  on  the  kth  row  of  particles.  For  a  straight 
crack,  U  and  G  are  given  by  equations  (2.1),  (2.3)  and  (2.14);  for  a  wavy  crack  they  may 
be  found  using  the  perturbation  procedure  outlined  in  the  previous  section.  For  particles 
perfectly  pinned  to  the  matrix,  ^fx^)  =  0.  For  debonding  particles,  S(x<rt' )  may  be  found 
by  eliminating  between  (2.49)  and  (2.51).  If  particles  at  x^^x^^.-.x^*  remain 
perfectly  pinned,  while  particles  at  x^N+1\ x^N^\ . .  .x^M^  are  debonding,  we  obtain  the 
following  expressions  for  P^  and  6^ 

M 

-U(x(n))  =  P{k) G(x(*},x(n))  1  <  n  <  N,  (2.52) 

k=l 

M  0  \P  JH. 

2AC/(x(n))-2AP0  Y  G(x(*\x(n))  =  6(x(n>)-=^  Y  6(x<*))G(x(*\x(n)) 

*=V+i  *=*+! 

N 

-2A^P(*)G(x(*),x("))  N<n<M.  (2.53) 

k=i 

Equations  (2.52)  and  (2.53)  may  be  solved  for  the  unknown  pW  over  the  perfectly  pinned 
grains,  and  6(x(fc))  at  debonding  particles.  P(t)  on  debonding  particles  is  then  deduced 
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using  (2.51).  Finally,  the  stress  intensity  factor  on  the  semi-infinite  crack  front  is  calculated 
using  (2.19). 

Although  equations  (2.52)  and  (2.53)  axe  linear,  the  number  of  perfecly  pinned  particles 
N  is  not  known  a  priori,  and  must  be  determined  using  an  iterative  procedure.  The 
following  technique  is  used  here.  Initially,  we  assume  *hat  all  particles  are  perfecly  pinned, 
and  solve  (2.52)  for  P^k\  The  condition  that  P^k^  <  Pq  is  then  checked:  any  grains  where 
this  criterion  is  violated  are  allowed  to  debond.  A  new  distribution  of  P  and  6  is  then 
calculated,  and  this  procedure  is  repeated  until  convergence  occurs.  Problems  involving 
frictional  behaviour  are  load-history  dependent.  The  solution  calculated  here  implicitly 
assumes  that  the  crack  is  subjected  to  a  monotonically  increasing  load,  so  that  the  force 
on  each  pinning  particle  never  reverses  sign. 

For  a  straight  crack,  the  influence  functions  G  and  H  are  known  in  closed  form.  For  a 
wavy  crack,  they  may  be  calculated  using  the  perturbation  procedure  outlined  in  section  2.4. 
However,  it  is  found  that  in  some  cases  over  500  rows  of  debonding  grains  may  be  left  in 
the  wake  of  the  crack:  a  prohibitive  amount  of  computing  time  is  required  to  update  the 
influence  functions  for  this  number  of  rows.  We  have  avoided  this  difficulty  by  updating  G 
and  H  explicitly  for  the  10  rows  of  particles  closest  to  the  crack  front,  and  calculating  the 
remainder  approximately.  The  function  G(x^n\x^^)  is  estimated  by  replacing  the  wavy 
crack  front  with  an  equivalent  straight  crack.  The  equivalent  crack  is  chosen  by  matching 
the  average  distance  between  the  row  of  particles  and  the  crack  front. 

Estimating  H  is  less  straightforward.  However,  inspection  of  equation  (2.38)  shows 
that  if  x  —  z(s)  >>  Sa(s),  the  incremental  expression  for  H  may  be  approximated  by 

1  fc/2 

6H(s,x)=—  /  H(A,x)D(A,s)[6a(A)  —  6a(s)n(X)  ■  n(s)l  dX.  (2.54) 

^  J  —  c/2 

Comparison  with  (2.37)  shows  that  this  is  equivalent  to  the  expression  for  6K(s),  with 
K( A)  =  H(A,x)  and  6K°°  =  0.  Let  Ho(A,x)  be  the  magnitude  of  H  for  a  straight  crack, 
as  defined  in  (2.18).  Noting  that  Ho(A,x)  is  almost  independent  of  A,  and  that  A'(A)  is 
constant,  we  conclude  by  induction  that  if  K{\)  is  known  for  the  wavy  crack  under  a  remote 
load  K°°,  H  may  be  deduced  from 

H(A,x)«H0(A,x)^.  (2.55) 

Essentially,  (2.55)  states  that  the  stress  intensities  induced  by  a  row  of  point  loads  applied 
a  long  way  from  the  front  of  a  wavy  crack  may  be  estimated  by  replacing  the  point  loads 
by  an  equivalent  uniform  remote  stress.  The  uniform  stress  is  chosen  to  induce  the  same 
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stress  intensity  factor  on  a  straight  crack  as  the  point  loads.  The  same  result  is  obtained 
by  matching  energy  release  rates  for  the  two  cases.  The  accuracy  of  the  estimates  of  G  and 
H  has  been  tested  by  comparing  the  approximate  values  with  the  exact  results  calculated 
using  the  finite  perturbation  method.  For  the  tenth  row  of  particles  behind  the  crack  front, 
the  error  in  the  approximate  solution  is  less  them  1%.  Any  arbitrary  number  of  debonding 
particles  may  thus  be  added  to  the  wake  of  the  crack. 

2.5)  Numericed  Implementation 

The  numerical  procedure  used  to  add  pinning  particles  to  the  crack  is  straightforward, 
and  merely  involves  solving  a  set  of  linear  equations.  However,  special  procedures  are  re¬ 
quired  to  solve  the  system  of  singular  integral  equations  that  arises  during  the  perturbation 
of  the  '"•ack  geometry.  The  numerical  technique  used  to  discretise  the  integral  equations 
follows  closely  the  procedure  given  in  Bower  and  Ortiz  (1990),  and  will  only  be  briefly 
summarised  here.  One  wavelength  of  the  crack  front  is  divided  into  m  elements,  with  three 
nodes  s2k-2,  $2k-i  and  s2k  on  the  &th  element.  In  addition,  2m—  1  collocation  points  such 
that  Sj  <  tj  <  Sj+ 1  are  defined.  The  crack  advance  is  then  approximated  using  Hermitian 
interpolation  between  values  of  Sa(s2k)  and  d6a(s2k)/ds  for  k  =  1,2. .  .m.  The  variation 
of  stress  intensity  factor  over  each  element  is  approximated  by  means  of  a  Lagrangian  in¬ 
terpolation  between  nodal  values;  similarly,  to  approximate  the  weight  functions  H,  D,  H 
and  D  we  write: 


H<-  ,  kc  -  ^5  +  A:c’x(n)) 

1  ^  ;  (z(s)  —  -1-  kLj]  •  [z(s)  —  x(n)  +  kLj]  ’ 

(2.56) 

n/  i  /  ^  d(s  -f  me,  t) 

(5  +  ,<}  ~  [z(s)  -  z(t)  +  kLj]  ■  [z(s)  -  z (t)  4  *Ij]’ 

(2.57) 

Hf„  (nK  _  h(s,x(n>)(sinh27T0)/(2-0) 

(T/7r)2(sin2  7rC  4  sinh2  7r^)  ’ 

(2.58) 

D/_  _  d(s,t)(sinh27r7)/(27T7) 

(L/7r)2(sin2  7T77  4  sinh2  717)  ’ 

(2.59) 

6  =  (xi  -  zi(s))/L;  (  =  (x2  ~z2(s))/L 

7  =  (*i(5)  -  zi (t))/L;  7]  =  (z2(s)  -  z2(t))/L 

and  use  three  point  Lagrangian  interpolation  between  nodal  values  to  approximate  the 
functions  h,  d ,  h  and  d.  The  functions  h(s  -f  kc,x)  and  d(s  4  kc,t )  need  to  be  stored  (and 
updated)  over  a  number  of  wavelengths  of  the  crack  front,  k  —  —M,  —M  4  1 ...  M  —  1,  M, 
as  discussed  in  Bower  and  Ortiz  (1990).  In  addition,  it  is  convenient  to  define  two  more 
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functions  /i00('S,x)  and  doo(s,<),  such  that 

hoo(s,x)  =  iim  h(s  ±  kc,  x),  (2.60) 

k  —  oc 

doc(s,t)  =  lim  d(s  ±  kc,  t).  (2.61) 

k— »oo 

The  infinite  series  in  equations  (2.18),  (2.39)  and  (2.40)  axe  then  evaluated  by  splitting 
each  series  into  a  periodic  and  a  non-periodic  part.  The  periodic  part  may  be  summed 
explicitly,  and  the  remaining  series  is  summed  numerically  by  truncating  it  at  a  finite 
number  of  terms.  For  example,  to  evaluate  the  series  in  (2.39),  we  write: 


d(s,t)(sinh27T7)/(27T7)  y>  d(s  +  kc,t) 

(L/n  )2(sin2  mj  +  sinh2  ny)  i-^oo  —  Xl(0]2  +  [^(-s)  —  ^2(0  +  kL}2 

_  y?  doo(s,t) _ y?  _ d(s  +  kc,  t)  -  doo(s,t) _ 

[ari(s)  -  Xi(t)]2  +  [x2(s)  -  x2(<)  +  kL}2  +  [xi(a)  -  xi(t)]2  +  [aj2(s)  -  x2(t)  +  kL}2 

d00(s,f)(sinh27r7)/(27r7)  ^  y^  _ d(s  +  kc,t)  -  doojs^) _ 

(L/7r)2(sin2  ttt]  +  sinh2  ^7)  [xi(s)  —  Xi(i)}2  +  [x2(s)  —  x2(t) -\- kL}2 


The  choice  of  a  suitable  value  for  M  is  discussed  in  Bower  and  Ortiz  (1990):  normally, 
M  =  2  is  sufficient.  The  series  in  (2.18)  and  (2.40)  may  be  evaluated  in  a  similar  way,  but 
for  reasons  of  space  will  not  be  written  out  in  full  here. 


To  solve  the  integral  equations  (2.45)  and  (2.46),  the  integrands  are  expressed  in  terms 
of  nodal  values  of  8a,  I\ ,  H  and  D  by  means  of  the  adopted  interpolation  functions.  As  a 
result,  the  integrals  may  be  evaluated  explicitly,  and  the  equations  may  be  written  in  the 
linear  form 

m 

0  =  J2  Wi(x^)6a(s2l)  +  lF,(x(n))<5a'(52t) 

«=o 

+  J^( Sp{k)  -  G(x(n)>x(<:))>  (2-63) 

k=  1  V  ' 

CT/CO  m 

=  lr=-*(‘>>  +  E  Vi(tj)Sa(s2i)  +  Vi(tj)6a'(s2i ) 

1=0 

+  Y,  (spii)  -  6-jS-p(k))  H(x(">,x(‘>),  (2.64) 

*=1  '  ' 


where  lTT,(x),  B^(x),  V,(t)  and  Vi(t)  are  integration  weights  depending  on  the  current 
crack  configuration.  The  procedure  used  to  calculate  these  functions  is  given  in  Bower  and 
Ortiz  (1990).  Equations  (2.63)  and  (2.64)  may  then  be  solved  for  the  nodal  values  of  8a 
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and  6a' .  Subsequently,  the  numerical  approximations  to  the  weight  functions  need  to  be 
updated:  this  is  simply  a  matter  of  evaluating  the  appropriate  principal  value  integrals  in 
(2.32)  and  (2.33):  an  appropriate  method  is  outlined  in  Bower  and  Ortiz  (1990). 

3.  Results  and  discussion 

The  perturbation  method  outlined  in  the  preceding  sections  may  be  used  to  analyse  the 
behaviour  of  a  crack  as  it  propagates  through  a  reinforced  brittle  solid.  Here,  we  consider 
an  initially  straight  semi-infinite  crack  propagating  through  a  brittle  matrix  material  of 
toughness  K “at.  The  matrix  is  reinforced  by  a  rectangular  distribution  of  tough  circular 
particles,  having  a  spacing  L,  radius  R  and  toughness  K^T .  The  ratio  of  particle  radius  to 
spacing  is  restricted  to  the  range  0  <  R/L  <  0.25,  to  be  consistent  with  the  assumptions 
made  in  Section  2.  The  solid  is  subjected  to  remote  loads,  acting  so  as  to  induce  a  mode  I 
stress  intensity  factor  on  the  crack.  The  magnitude  of  the  load  is  parameterised  by  a 
remotely  applied  stress  intensity  factor  K°°,  as  defined  in  equation  (2.26).  Initially,  we 
assume  that  the  particles  are  perfectly  bonded  to  the  matrix  material,  and  fracture  in  the 
wake  of  the  crack.  The  analysis  will  subsequently  be  extended  to  include  the  effects  of 
debonding  between  particles  and  the  matrix. 

3.1  Perfectly  bonded  particles 

( a )  Crack  Trapping 

Consider  first  the  behaviour  of  a  crack  with  no  pinning  particles  in  its  wake,  as  it  meets  a 
row  of  tough  particles.  The  crack  may  surmount  the  row  of  obstacles  in  one  of  two  ways.  If 
the  toughness  of  the  particles  exceeds  that  of  the  matrix  by  only  a  small  amount,  they  are 
penetrated  by  the  crack,  which  propagates  through  them.  Alternatively,  if  the  K£*T  / 
exceeds  a  critical  value  (to  be  determined),  the  segments  of  crack  front  bowing  out  between 
the  particles  coalesce  with  one  another,  and  the  particles  are  left  intact  in  the  wake  of  the 
crack. 

The  latter  behaviour  is  illustrated  on  Figs  4  to  6.  Fig.  4  shows  the  shape  of  one 
wavelength  of  the  crack  front  as  it  propagates  around  the  first  row  of  tough  particles,  for 
a  particle  spacing  L/R  =  10.  The  front  is  initially  straight,  and  the  wake  of  the  crack 
is  free  of  pinning  particles.  When  the  remote  load  reaches  a  magnitude  K00  —  A™at , 
the  crack  starts  to  propagate  through  the  solid.  Parts  of  the  crack  front  contacting  the 
particles  arrest,  while  the  remainder  of  the  front  bows  out  between  the  particles.  Its  shape 
is  determined  by  the  condition  that  K(s)  =  A'“at  over  propagating  regions,  while  6a(s)  =  0 
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over  regions  in  contact  with  the  tougher  material.  Due  to  the  change  in  crack  geometry, 
the  stress  intensity  factor  increases  over  regions  of  the  crack  which  contact  the  particles: 
the  resulting  distribution  of  K(s )  along  the  crack  front  is  shown  in  Fig.  5,  at  various  stages 
of  growth.  The  results  are  shown  as  a  function  of  s/c,  where  the  arc-length  s  is  measured 
from  the  point  of  first  contact  of  the  crack  front  with  the  obstacle,  and  c  is  the  arc-length 
of  one  period. 

The  variation  of  remote  stress  as  the  crack  propagates  around  the  row  of  particles  is 
shown  in  Fig.  6.  The  load  is  plotted  as  a  function  of  (a) / R:  the  parameter  (a)  is  a  measure 
of  the  average  crack  advance,  defined  as  Af L,  where  A  is  the  area  swept  by  one  wavelength 
of  the  advancing  crack.  The  load  reaches  a  maximum  value  of  1.68,  and  subsequently 
decreases  as  areas  of  the  crack  front  ahead  of  the  particles  start  to  attract  one  another.  As 
the  load  decreases,  parts  of  the  crack  front  arrest,  and  the  propagating  region  progressively 
diminishes  in  size  to  a  small  zone  around  the  point  of  coalescence.  The  finite  perturbation 
method  cannot  be  used  to  analyse  in  detail  the  coalescence  of  the  crack,  due  to  the  singular 
behaviour  of  the  influence  function  D(s,t),  as  discussed  in  section  2.  The  calculation  was 
stopped  when  the  crack  front  reached  the  last  profile  shown  on  Fig.  4.  However,  the  remote 
load  must  fall  to  zero  at  the  point  of  coalescence:  since  the  stress  intensity  factor  is  singular 
where  the  crack  front  touches  itself,  it  follows  that  K°°  — ►  0  to  satisfy  the  condition  that 
K(s)  —  K ™at  over  the  propagating  region. 

In  practice,  the  coalescence  of  the  crack  can  only  be  a  quasi-static  process  as  shown  in 
Figs  4  and  6  if  the  solid  is  subjected  to  displacement  boundary  conditions  at  infinity.  Under 
a  steadily  increasing  remote  stress,  the  configuration  of  the  crack  will  become  unstable  when 
the  load  reaches  the  peak  value  I\°° /A'™at  =  1.68,  and  the  crack  will  propagate  dynamically 
until  it  is  arrested  by  the  next  row  of  particles.  This  behaviour  is  known  as  ‘pop-in’,  and 
has  been  observed  during  tests  on  Alumina  by  Swanson  et  al  (1987). 

Fig.  5  shows  that  the  stress  intensity  on  the  region  of  the  crack  in  contact  with  the  parti¬ 
cle  reaches  a  maximum  value  K(s )  =  2.84  at  s  =  0.  Consequently,  if  Kf ar/ A'™al  >  2.S4,  the 
crack  bypasses  the  first  row  of  obstacles,  leaving  a  row  of  pinning  particles:  its  subsequent 
propagation  is  discussed  in  detail  in  the  next  section.  Alternatively,  if  A'?ar/A'™at  <  2.84, 
the  crack  front  penetrates  the  particles.  In  this  case,  the  obstacles  fracture,  and  no  bridging 
particles  are  formed.  The  behaviour  of  a  crack  breaking  through  a  row  of  obstacles  has 
been  discussed  in  detail  by  Gao  and  Rice  (1990).  They  distinguish  between  two  different 
cases:  stable  (or  regular),  and  unstable  (irregular)  penetration.  Stable  breakthrough  im¬ 
plies  that  the  crack  fully  penetrates  the  obstacles  before  the  maximum  load  is  reached. 
The  penetration  is  unstable  if  the  peak  load  occurs  when  the  particles  are  only  partially 
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penetrated. 

The  maximum  load  under  conditions  of  stable  penetration  occurs  when  the  particles  are 
fully  penetrated,  and  the  line  fraction  of  the  crack  front  inside  the  obstacles  is  a  maximum. 
In  this  case,  the  toughening  due  to  trapping  can  be  calculated  explicitly,  by  means  of 
equation  (2.26)  (Rose,  1975).  Since  K(x 2)  =  Kf™  for  \i2\  <  R  and  K(x 2)  =  A'™at 
otherwise,  we  immediately  obtain 

j  j2 

A' 00  =  j(Acmat)2  +  ^  ((Kcpar)2  ~  (A“al)2)  j  .  (3.1) 

This  expression  is  exact,  and  is  valid  provided  the  crack  can  reach  the  limiting  configuration. 
It  gives  an  upper  bound  to  the  toughening  due  to  trapping  if  the  crack  bypasses  the 
obstacles,  or  in  the  event  of  unstable  penetration. 

The  shape  of  the  crack  at  the  maximum  load  may  be  found  using  the  finite  perturbation 
procedure.  The  required  configuration  is  determined  by  the  condition  that  K(s)  =  A'Par 
inside  the  particles,  and  K(s)  =  A'™at  outside.  This  can  be  satisfied  by  applying  a  succes¬ 
sion  of  perturbations  such  that  6K(s)  =  constant  for  1 0*2 1  <  A,  and  6I\(s )  ~  0  otherwise. 
The  change  in  crack  geometry  may  be  found  by  solving  equation  (2.31)  for  Sa(s)  at  each 
increment.  In  this  wav,  the  limiting  crack  configuration  may  be  found  for  a  range  of  val¬ 
ues  of  A'Par  without  having  to  solve  explicitly  for  the  load  history  in  each  case.  A  set 
of  results  is  shown  in  Fig.  7,  for  a  particle  spacing  L/R  =  6.67.  It  can  be  seen  that  if 
A'Par/A’^al  >  1.7,  the  limiting  crack  shape  cannot  be  contained  within  the  particles.  Un¬ 
der  these  conditions,  the  maximum  load  occurs  before  the  obstacle  is  fully  penetrated.  This 
case  has  been  classified  as  unstable  penetration  by  Gao  and  Rice  (1990). 

If  the  particles  are  not  fully  penetrated  at  maximum  load,  equation  (3.1)  does  not 
apply,  and  the  toughness  of  the  reinforced  solid  must  be  calculated  numerically.  A  set  of 
crack  profiles  for  an  unstably  penetrated  obstacle  is  shown  on  Fig.  8,  for  A’Par/A'f al  =  2.4 
and  L/R  =  6.667.  The  maximum  load  was  found  to  be  A' 00 / A'™at  =  1.53:  the  correspond¬ 
ing  value  predicted  by  equation  (3.1  )is  1.55  Tu?  crack  front  bowing  between  the  particles 
joins  up  on  itself  after  the  obstacles  are  penetrated  and  a  row  of  partially  penetrated  par¬ 
ticles  is  left  in  the  wake  of  the  crack.  Its  subsequent  behaviour  is  discussed  below.  For 
small  panicles  ( R/L  <  0.18),  the  penetrated  obstacles  fracture  before  the  crack  reaches  the 
second  row  of  particles.  Under  these  conditions  crack  bridging  does  not  contribute  to  the 
toughness  of  the  solid:  the  maximum  load  occurs  while  the  crack  penetrates  a  single  row 
of  obstacles,  at  the  point  of  instability.  However,  if  R/L  >  0.18,  the  penetrated  particles 
remain  intact.  This  situation  cannot  be  modelled  exactly  using  the  present  method,  since 
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the  pinned  regions  in  the  crack  wake  have  been  assumed  to  be  circular.  The  penetrated  par¬ 
ticles  are  therefore  replaced  by  a  row  of  equivalent  circular  obstacles,  chosen  by  matching 
areas. 

These  calculations  have  been  repeated  for  a  range  of  values  of  L/R  and  A^Par/A'™al. 
The  results  are  summarised  on  Figs.  9  and  10.  Fig.  9  shows  the  mechanism  of  trapping  as  a 
function  of  the  particle  toughness  and  spacing.  In  the  regime  marked  ‘stable  penetration,’ 
the  crack  front  fuliy  penetrates  the  obstacles  at  the  maximum  load,  and  the  toughness 
is  given  by  equation  (3.1).  For  slightly  larger  values  of  AfPar / A'“at  in  the  ‘unstable  pen¬ 
etration’  region,  the  maximum  load  occurs  before  the  obstacles  are  fully  penetrated.  In 
the  region  marked  ‘unstable  bridging’,  a  row  of  partially  penetrated  particles  remains  in 
the  wake,  but  fractures  before  the  crack  front  reaches  the  second  row  of  obstacles.  In 
the  ‘Partial  Bridging’  regime,  the  penetrated  particles  remain  intact  and  contribute  to  the 
toughness  of  the  solid  by  pinning  the  crack  faces.  Finally,  if  Kf**  /  K™'x  lies  in  the  area 
marked  ‘particles  bypassed.’  the  obstacles  are  left  intact  in  the  wake  of  the  crack.  The 
subsequent  behaviour  of  the  crack  is  discussed  in  more  detail  below.  Fig.  9  shows  that  the 
critical  value  of  A'|>ar//A'fiat  required  for  bridging  particles  to  be  formed  increases  as  L/R 
is  reduced.  If  there  is  a  large  distance  between  the  particles,  the  crack  tends  to  bow  out 
further,  so  that  the  crack  is  more  likely  to  coalesce  ahead  of  the  particles. 

The  effective  toughness  of  the  reinforced  material  is  shown  as  a  function  of  A'Par / K™3,1 
on  Fig.  10.  for  various  values  of  L/R.  The  toughness  enhancement  due  to  crack  trapping 
is  a  monotonically  increasing  function  of  A'far/A'™at  and  R/L.  Even  if  no  crack  bridging 
takes  place,  the  particles  can  double  toughness  of  a  brittle  matrix  by  trapping  the  crack 
front. 

(b)  Bridging  combined  with  trapping 

If  the  particle  toughness  lies  in  the  regime  marked  ‘particles  bypassed’  in  Fig.  9.  at 
least  one  row  of  bridging  particles  will  be  formed  in  the  wake  of  the  crack.  This  situation 
is  analysed  by  adding  pinning  particles  behind  the  crack  front  and  subsequently  modifying 
its  shape,  as  discussed  in  section  2.2.  The  propagation  of  the  crack  immediately  after 
coalescence  cannot  be  analysed  using  this  procedure,  since  the  calculation  assumes  that 
the  distance  between  the  pinning  panicles  and  the  crack  front  greatly  exceeds  the  particle 
radius.  We  therefore  resume  the  analysis  after  the  semi-infinite  crack  has  propagated  a 
short  distance  beyond  the  first  row  of  pinning  particles.  Fig.  12  shows  the  load  as  a  function 
of  crack  advance  while  the  crack  bypasses  several  rows  of  particles:  the  dashed  portions  of 
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this  curve  occur  during  coalescence  and  have  not  been  calculated  explicitly.  We  have  made 
use  of  the  result  that  K°°  must  fall  to  zero  at  the  point  of  coalescence  in  order  to  deduce 
the  probable  shape  of  the  K°°  —  (a)  curve  in  this  regime.  Similarly,  the  dashed  regions  on 
Fig.  13  have  been  deduced  by  noting  that  ATpar  — ►  0  at  coalescence. 

A  set  of  crack  front  profiles  during  the  formation  of  four  pinning  particles  is  shown  on 
Fig.  11.  It  can  be  seen  that  while  the  crack  propagates  between  the  first  and  second  rows  of 
particles,  its  front  remains  approximately  straight.  At  this  stage,  the  remote  load  and  the 
stress  intensity  around  the  pinning  particles  axe  closely  approximated  by  the  results  for  a 
straight  crack  bridged  by  single  row  of  obstacles,  given  in  equations  (2.23)  and  (2.24).  It  is 
straightforward  to  show  that  K°°  / given  by  (2.23)  is  a  decreasing  function  of  the  crack 
advance  b  for  all  L/R.  Until  the  crack  front  is  trapped  by  the  second  row  of  particles,  the 
remote  load  therefore  always  decreases  as  the  crack  advances.  This  situation  is  similar  to 
a  crack  propagating  through  a  material  of  decreasing  toughness:  by  differentiating  (2.23), 
the  rate  of  change  of  effective  toughness  with  crack  advance  is  found  to  be 


dKf 
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— AT™at  cosh  zb/ L 
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where  b  is  the  distance  from  the  particles  to  the  crack  front,  and  F(b/L)  is  defined  in  (2.25). 
For  b/L  in  the  range  0.5  <  b/L  <  1.5,  this  may  be  approximated  by 
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Bower  and  Ortiz  (1990)  have  investigated  the  behaviour  of  a  straight  semi-infinite  crack 
propagating  through  material  of  decreasing  toughness.  They  have  shown  that  the  con¬ 
figuration  of  the  crack  is  not  unique,  in  that  the  crack  front  may  not  remain  straight.  In 
particular,  small  imperfections  in  the  geometry  of  the  crack  containing  a  critical  wavelength 
ccrjt  are  likely  to  grow  as  the  crack  advances.  The  critical  wavelength  is  given  by 
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For  the  case  in  hand 


b(L/R  -  0.93 L/bf 
6.5 


This  result  suggests  that  shortly  after  the  crack  bypasses  a  row  of  obstacles  ( b  small),  the 
crack  front  may  adopt  a  wavy  configuration,  and  advance  by  forming  a  series  of  short 
wavelength  ‘fingers’.  However,  for  the  range  of  L/R  considered  here,  ccrjt  >  L  for  b  ss  L. 
The  present  analysis  constrains  the  crack  front  to  remain  periodic  with  wavelength  L,  so 
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once  the  distance  from  the  bridging  particles  to  the  crack  front  is  comparable  to  the  obstacle 
spacing,  the  crack  front  remains  approximately  straight. 

The  stress  intensity  factor  on  the  bridging  particles  is  plotted  as  a  function  of  crack 
advance  on  Fig.  13,  for  L/R  =  10.  While  the  crack  propagates  between  the  first  and 
second  rows,  ATpar/AT“at  is  closely  approximated  by  the  result  for  a  straight  crack,  given 
in  equation  (2.23).  The  particles  fracture  when  ifpar  =  ATpar:  this  condition  may  be  used 
to  calculate  the  value  of  ATpar  required  for  the  obstacles  to  remain  intact  until  the  crack 
front  reaches  the  second  row  of  obstacles.  Putting  6  «  L  in  (2.23)  we  find  that: 

ifpar  _  rr  1 

tfmat  ~  V  2  R  7r/4  -  0.729 R/L '  ^3'6^ 

For  large  L/R,  the  particle  toughness  for  the  pinning  regions  to  be  stable  exceeds  the 
critical  value  required  for  the  crack  to  bypass  the  particles.  For  example,  if  L/R  =  10, 
the  crack  bypasses  the  first  row  of  particles  if  ATpar/AT™at  >  2.79.  Equation  (3.6)  shows 
that  the  particles  fracture  before  the  crack  front  reaches  the  second  row  of  obstacles  unless 
A'par /K™3t  >  3.14.  If  A'par / Ar™at  lies  between  2.79  and  3.14,  the  crack  bypasses  the  first 
row  of  obstacles,  which  subsequently  fractures  while  the  crack  advances  towards  the  second 
row. 

If  Apar/A’®at  exceeds  3.14,  the  bridging  particles  remain  intact  until  the  crack  is 
trapped  by  the  second  row  of  obstacles.  At  this  point,  the  remote  load  starts  to  increase 
while  the  crack  bows  out  between  the  obstacles,  as  shown  in  Fig.  12.  At  the  same  time,  the 
stress  intensity  on  the  bridging  particles  increases  rapidly,  Fig  13.  Eventually,  K°°  exceeds 
the  maximum  value  reached  while  bypassing  the  first  row  of  obstacles:  for  L/R  —  10,  the 
critical  value  A”00  =  1.43  is  reached  when  ( a)/R  =  11.  At  this  point  A’par/A'“at  =  3.54. 
Consequently,  if  the  particle  toughness  lies  in  the  range  2.84  <  A'par/A’fiat  <  3.54,  the 
maximum  load  occurs  at  the  instant  when  the  crack  bypasses  the  first  row  of  particles. 
Although  bridging  particles  are  subsequently  formed,  they  fracture  at  a  lower  load,  and  do 
not  contribute  to  the  toughness  of  the  reinforced  solid  This  behaviour  will  be  referred  to 
as  ‘unstable  bridging.’  The  range  of  Apar/A'™at  for  unstable  bridging  to  occur  is  shown  as 
a  function  of  R/L  on  Fig.  9.  In  this  regime,  variations  in  particle  toughness  produce  no 
change  in  the  effective  toughness  of  the  reinforced  material. 

A  similar  argument  applies  to  partially  penetrated  obstacles  in  the  wake  of  the  crack. 
If  the  obstacles  are  small,  they  fracture  before  the  crack  front  reaches  the  second  row  of 
particles.  Large  particles  remain  intact  for  long  enough  to  contribute  to  the  toughness  of 
the  reinforced  solid  by  pinning  the  crack  faces. 
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The  behaviour  of  the  crack  as  it  bypasses  further  rows  of  particles  is  similar.  As  the 
particle  toughness  is  increased,  several  stable  rows  of  pinning  particles  may  form  in  the  wake 
of  the  crack.  Fig.  9  shows  the  number  of  bridging  particles  as  a  function  of  K%*T  /  K™*1  and 
L/R.  In  the  regions  marked  ‘unstable  bridging,’  an  additional  row  of  particles  is  formed 
but  fractures  without  improving  the  toughness.  Under  these  conditions,  the  maximum  load 
is  reached  at  the  point  where  the  crack  just  bypasses  a  new  row  of  obstacles.  Otherwise, 
the  maximum  load  is  reached  when  the  row  of  bridging  particles  furthest  from  the  crack 
front  fractures.  In  the  steady  state,  bridging  particles  fracture  at  the  same  rate  in  the  wake 
of  the  crack  as  they  axe  formed  at  its  front.  The  load  fluctuates  considerably  as  the  crack 
propagates:  this  is  a  consequence  of  the  periodicity  of  the  distribution  of  particles.  Smaller 
variations  in  load  would  be  expected  if  the  particle  distribution  were  random. 

The  maximum  remote  load  has  been  plotted  as  a  function  of  the  particle  toughness 
on  Fig  14,  for  various  values  of  R/L:  K°°  may  be  interpreted  as  the  effective  toughness 
of  the  reinforced  solid.  It  can  be  seen  that  for  A£ar /AT“at  >  7,  K°° /K™*1  varies  linearly 
with  APar/A™at,  and  is  closely  approximated  by  the  relation: 

K°°  nR  A7ar 

— — -  =  3.09——^—.  (3.7) 

Kf al  L  A™at  v  ' 

It  is  not  difficult  to  explain  this  behaviour.  Consider  a  straight,  semi-infinite  crack  propa¬ 
gating  through  a  material  of  toughness  Kf1*1  which  contains  a  rectangular  distribution  of 
particles,  toughness  A’(Par,  radius  R  and  spacing  L.  As  the  crack  propagates  a  distance  L, 
the  work  done  per  unit  area  in  fracturing  the  particles  varies  as  7rA2(A'|>ar)2/X2,  while  the 
work  done  in  cleaving  the  matrix  is  proportional  to  (1  —  7rA2/X2)(A'Jnat)2.  It  follows  that 
the  effective  toughness  of  the  composite  material  is  given  by 


Kf  I  R2  /  (A?ar)2 

ATcmat  ~  I/  +?rX2  V(A'cmat)2 


Expanding  (3.8)  for  K?*T  j  K™3-*  >>  1  shows  that 


A™at 


r-R  Kr 

V  '  l  A^mat’ 


giving  the  linear  behaviour  observed  in  Fig.  14.  The  constant  coefficients  in  (3.7)  and 
(3.9)  differ,  since  equation  (3.9)  gives  the  average  effective  toughness,  while  (3.7)  gives  the 
maximum  value,  and  includes  the  transient  effect  of  crack  trapping. 

It  is  evident  that  adding  small  quantities  of  tough  particles  to  a  brittle  material  earn 
substantially  improve  its  toughness.  While  the  toughness  enhancement  due  to  crack  tip 
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microcracking  or  deflection  is  modest  at  best,  the  present  analysis  shows  that  crack  trapping 
combined  with  bridging  can  improve  the  fracture  toughness  by  a  factor  of  5  or  more. 
However,  several  criteria  must  be  satisfied  in  order  to  achieve  these  levels  of  toughness. 
The  particles  must  be  perfectly  bonded  to  the  matrix  material,  and  their  toughness  must 
be  of  the  order  of  10  times  that  of  the  matrix.  This  may  be  the  case  in  ceramics  and  glasses 
reinforced  by  metallic  particles,  but  is  unlikely  in  monolithic  ceramics.  The  analysis  also 
shows  that  only  a  very  small  number  of  perfectly  pinned  particles  are  likely  to  form  in  the 
wake  of  a  crack.  Even  if  ivPar/i^“at  ss  15,  only  10  particles  remain  intact  for  L/R  =  20. 

(c)  Comparison  with  experiments 


Krstic  et  al  (1981)  have  measured  the  toughness  of  glass  reinforced  by  partially  oxidized 
aluminium  particles.  The  assumptions  made  in  the  present  analysis  satisfy  the  conditions 
of  their  experiments.  The  oxide  layer  on  the  particles  forms  a  strong  bond  with  the  AI2O3 
in  the  glass,  and  the  elastic  constants  and  thermal  expansion  coefficient  of  both  materials 
are  closely  matched.  The  measured  toughness  (normalised  by  Kjc  for  pure  glass)  is  shown 
on  Fig.  14,  as  a  function  of  the  volume  fraction  of  particles,  V/. 


In  order  to  compare  experimental  and  theoretical  results,  it  is  necessary  to  relate 
the  particle  spacing  used  in  the  analysis  to  an  equivalent  volume  fraction  of  randomly 
distributed  particles.  It  is  straightforward  to  show  that  if  a  plane  cuts  through  a  random 
distribution  of  spheres,  the  area  fraction  of  the  plane  inside  the  particles  is  equal  to  their 
volume  fraction.  Matching  area  fractions  for  the  periodic  and  random  distributions,  we 


find  that 


(3.10) 


Using  (3.10),  equation  (3.7)  may  be  expressed  in  terms  of  Vf  as 


K°° 

A™at 


- h'Par 


(3.11) 


This  result  has  been  plotted  on  Fig.  14.  Taking  A'p^/A'™1*  =11  gives  excellent  agreement 
between  theory  and  experiment.  However,  this  choice  of  particle  toughness  is  arbitrary,  and 
is  probably  an  underestimate.  With  higher  values  of  A'Par,  equation  (3.11)  overestimates  the 
toughness  of  the  composite:  this  is  to  be  expected,  since  the  periodic  particle  distribution 
used  in  the  theory  overestimates  the  toughness  enhancement  due  to  crack  trapping  by  a 
random  distribution. 
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The  R-curve  behaviour  of  a  monolithic  ceramic  has  been  measured  by  Vekinis  et 
al  (1990).  Their  results  axe  shown  on  Fig.  16.  The  toughness  reaches  a  steady  state 
after  11mm  of  crack  growth,  suggesting  that  the  length  of  the  bridged  zone  exceeds  2000 
grain  diameters.  According  to  the  present  analysis,  particles  would  need  to  be  over  100 
times  tougher  than  the  matrix  to  remain  intact  over  this  length:  equation  (3.7)  would 
then  predict  a  corresponding  toughness  enhancement  of  the  order  of  a  factor  of  30.  It 
is  clearly  not  reasonable  to  assume  that  grains  remain  perfectly  bonded  to  the  matrix  in 
monolithic  ceramics.  The  R-curve  behaviour  observed  in  Fig.  16  is  more  likely  to  be  caused 
by  frictional  crack  bridging.  This  issue  is  addressed  in  the  next  section. 

3.2  Frictional  crack  bridging 

If  particles  axe  weakly  bonded  to  the  matrix,  the  interface  between  particles  and  matrix 
fails  before  the  particles  fracture,  so  that  they  are  pulled  out  in  the  wake  of  the  crack.  The 
idealised  model  of  this  process  was  discussed  in  section  2.4.  We  assume  each  grain  to  be 
cylindrical,  having  a  radius  R  and  height  6crit,  with  its  axis  perpendicular  to  the  plane  of 
the  crack.  The  particles  are  assumed  to  be  sufficiently  tough  so  as  not  to  be  penetrated  by 
the  crack.  The  foremost  row  of  particles  remains  bonded  to  the  matrix,  so  that  the  effects 
of  crack  trapping  are  as  discussed  in  the  preceding  sections:  the  particles  axe  bypassed  by 
the  crack  front  and  remain  intact  in  the  wake  of  the  crack.  The  particles  remain  perfectly 
bonded  until  the  pinning  force  on  a  row  of  grains  exceeds  P0  =  ^RScritf^r,  at  which  point 
they  start  to  be  pulled  free  from  the  matrix.  Thereafter,  the  pinning-force  is  related  to  the 
crack  opening  displacement  by  equation  (2.49). 

The  toughening  due  to  frictional  crack  bridging  can  be  expressed  in  terms  of  four 
non-dimensional  groups 

K°°  _  (  PoScrit  R  L26cth  6ciit\ 

K^t  r  \L2A(K™')2'  V  ARP0  ’  R  )  '  1 ’ 

Here,  the  first  group  is  the  ratio  of  the  frictional  energy  dissipation  rate  to  the  critical  energy 
release  rate  of  the  matrix  material.  This  is  equivalent  to  the  parameter  (Kf&T / A”™at)2  in 
the  case  of  perfectly  pinned  particles.  The  second  group,  L/R  is  a  measure  of  the  volume 
fraction  of  particles,  as  before.  The  third  group  is  new,  and  requires  some  explanation.  The 
ratio  Po /L28CT\t  is  a  measure  of  the  slope  of  the  P  —  6  curve  during  pull-out  of  the  particles. 
In  contrast,  A R  is  related  to  the  compliance  of  the  faces  of  the  crack  under  a  distribution 
of  pressure.  Thus,  the  third  dimensionless  group  in  (3.12)  is  a  measure  of  the  ratio  of 
the  compliance  of  the  frictional  pinning  force-pull  out  curve  to  the  elastic  compliance  of 
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the  crack  faces.  This  parameter  has  an  important  influence  on  the  length  of  the  bridging 
zone.  The  last  dimensionless  group  governs  the  three-dimensional  shape  of  the  particles: 
for  simplicity,  we  will  assume  that  <5crit/i2  =  2  for  all  calculations  here. 

Typical  values  for  the  parameters  in  this  model  are  given  in  Table  1,  chosen  to  be 
representative  of  alumina.  Parameters  such  as  the  grain  size  R  and  the  elastic  constants  are 
easily  determined,  but  the  parameters  governing  frictional-pull  out  can  only  be  estimated. 
Here,  we  have  assumed  that  the  compressive  residual  stress  pinning  each  grain  to  the  matrix 
or  &  120MPa.  Stresses  of  this  magnitude  have  been  predicted  by  Ortiz  and  Suresh  (1989) 
and  measured  by  Blendel  and  Coble  (1982).  The  coefficient  of  friction  has  been  taken  as  0.6. 
For  a  grain  size  of  the  order  of  4/zm,  this  leads  to  a  maximum  pull-out  force  Po  ~  0.008N. 

Fig.  15  shows  the  predicted  R-curve  behaviour  due  to  frictional  crack  bridging,  for 
an  appropriate  value  of  Potent / L2 A(JsTf'at)2 .  The  fluctuation  in  load  as  the  crack  bypasses 
individual  rows  of  particles  has  not  been  shown  in  Fig  15:  we  have  assumed  that  the  solid  is 
subjected  to  a  monotonically  increasing  remote  stress,  and  ignored  unstable  crack  growth. 
The  effective  toughness  was  calculated  by  adding  rows  of  particles  to  the  wake  of  the  crack, 
and  then  finding  the  maximum  value  of  remote  load  as  the  crack  bypasses  the  next  row. 
This  procedure  gives  an  upper  bound  to  the  combined  toughening  due  to  bridging  and 
trapping. 

The  R-curve  initially  rises  steeply  as  the  crack  bypasses  the  first  row  of  particles. 
The  crack  growth  during  the  initial  increase  in  stress  is  of  the  order  of  the  grain  size,  so 
it  is  unlikely  that  this  part  of  the  curve  could  be  detected  experimentally.  Instead,  the 
initial  toughness  of  the  material  would  appear  to  be  greater  than  that  of  the  matrix.  The 
initial  toughness  of  alumina  measured  by  Vekinis  et  al  (1990)  (see  Fig.  15)  is  of  the  order 
of  3.3  MPay'm:  the  inherent  matrix  toughness  to  give  this  initial  value  is  l^MPa-^/m. 
The  equivalent  surface  energy  of  6  Jm-2  is  in  good  agreement  with  measured  values  for 
Alumina.  Crack  trapping  can  therefore  account  for  the  observed  discrepancy  between  the 
surface  energy  and  fracture  toughness  of  polycrystalline  alumina.  However,  the  origin  of 
the  tough  grains  responsible  for  trapping  is  yet  to  be  determined  with  certainty. 

Once  the  particles  have  been  bypassed,  they  remain  perfectly  bonded  to  the  matrix 
for  only  a  short  time.  For  the  parameters  used  here,  no  more  than  a  single  row  was  found 
to  remain  fully  pinned  The  remaining  particles  pull  out  in  the  wake  of  the  crack.  As  a 
result,  the  toughness  enhancement  due  to  bridging  is  greatly  reduced:  although  over  500 
particles  remain  in  the  bridged  zone,  the  toughness  enhancement  is  only  of  the  order  of 
20-30%.  The  predicted  increase  in  toughness  is  in  good  agreement  with  the  experimentally 
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measured  values  shown  in  Fig.  15;  furthermore,  the  measured  and  calculated  length  of  the 
bridging  zones  are  comparable. 

The  effects  of  varying  the  parameters  P0(5crit/L2A(i;fJnat)2,  R/L  and  L26Cjh/ARPq  are 
shown  in  Figs.  16  to  18.  Fig.  16  shows  the  influence  of  the  frictional  energy  dissipation 
on  the  toughness  increase  due  to  bridging  A Kr,  while  Fig.  17  shows  the  variation  of  the 
length  of  the  bridged  zone.  The  compliance  ratio  L28CTlt/ARPo  has  little  influence  on  the 
toughness:  varying  this  parameter  from  10  x  103  to  100  x  103  produced  only  a  5%  change  in 
toughness.  However,  the  length  of  the  bridged  zone  is  controlled  by  the  compliance  ratio, 
as  seen  on  Fig.  18.  This  suggests  that  the  toughness  is  mainly  determined  by  the  area 
under  the  pinning  force-pull  out  curve,  and  is  not  greatly  affected  by  its  shape. 

The  influence  of  grain  size  on  the  toughness  of  Alumina  is  an  intriguing  feature  of  the 
experimental  study  by  Vekinis  et  al.  (1990).  The  toughness  was  observed  to  increase  with 
grain  size,  but  the  length  of  the  bridging  zone  decreased.  The  effects  of  variations  in  grain 
size  may  be  calculated  using  the  theoretical  model  developed  here.  From  equation  (2.50), 
we  find  that  Pq  varies  as  R2.  If  the  volume  fraction  of  tough  grains  (parameterised  by 
L/R)  remains  fixed,  and  the  aspect  ratio  8Cnt/R  is  unchanged,  then  the  first  dimensionless 
group  in  (3.12)  varies  linearly  with  grain  size  R ,  while  the  other  groups  are  unaffected. 
The  effect  of  grain  size  on  toughness  and  bridging  length  may  therefore  be  deduced  from 
Figs.  16  and  17.  The  measured  toughness  agrees  well  with  the  calculated  values,  but  the 
theory  predicts  that  the  bridging  length  increases  with  grain  size,  contrary  to  experimental 
observations. 

The  reason  for  this  discrepancy  is  not  clear.  It  is  not  difficult  to  adjust  the  parameters 
in  the  theoretical  model  to  produce  a  simultaneous  increase  in  toughness  and  decrease  in 
bridging  length:  fox  example,  the  compliance  ratio  or  the  aspect  ratio  of  the  grains  could  be 
chosen  appropriately.  The  annealing  process  used  to  increase  the  grain  size  may  alter  the 
residual  stress  distribution  in  the  solid,  or  change  the  interfacial  toughness.  This  influences 
the  number  of  grains  left  intact  in  the  wake  of  the  crack.  However,  the  volume  fraction  and 
aspect  ratio  of  tough  grains  are  difficult  to  measure,  and  no  data  are  available  at  present 
to  allow  a  complete  study. 

4.  Conclusions 

A  three-dimensional  model  has  been  developed  to  calculate  the  toughness  of  a  brittle 
material  reinforced  by  tough  particles.  The  reinforced  material  has  been  idealised  as  an 
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elastically  homogeneous  solid,  containing  a  regular  distribution  of  tough  circular  particles. 
The  analysis  predicts  the  behaviour  of  an  initially  straight,  planar,  semi-infinite  crack  as 
it  propagates  through  the  solid  under  a  uniform  remote  stress.  The  particles  may  inhibit 
fracture  by  three  mechanisms:  firstly,  the  particles  trap  the  crack  front,  causing  it  to  bow 
out  between  them.  Secondly,  the  particles  may  remain  intact  in  the  wake  of  the  crack, 
thereby  pinning  its  faces  and  carrying  part  of  the  load.  Finally,  the  particles  may  debond 
from  the  matrix  and  dissipate  energy  in  frictional  sliding  as  they  are  pulled  out  in  the  wake 
of  the  crack.  All  three  mechanisms  have  been  investigated. 

The  results  for  particles  perfectly  bonded  to  the  matrix  are  summarised  in  Figs.  9 
and  10.  Fig.  9  shows  the  number  of  pinning  particles  that  form  in  the  wake  of  the  crack, 
as  a  function  of  the  obstacle  spacing  and  toughness.  If  exceeds  X“at  by  only  a  small 
amount,  the  obstacles  are  penetrated  by  the  crack  front,  and  no  bridging  particles  are 
formed.  Crack  trapping  is  then  the  only  toughening  mechanism.  If  the  particle  toughness 
exceeds  a  critical  value,  they  remain  intact  in  the  wake  of  the  crack.  The  toughness  en¬ 
hancement  due  to  combined  crack  trapping  and  bridging  is  shown  in  Fig.  10.  Asymptotic 
expressions  have  been  found  for  the  toughness  of  the  reinforced  solid,  both  for  particles 
whose  toughness  is  comparable  to  that  of  the  matrix,  and  also  for  particles  much  stronger 
than  the  matrix.  If  A”Par  <  3 AT™at,  the  obstacles  are  penetrated  by  the  crack  front,  and 
no  bridging  particles  form  in  the  crack  wake.  Under  these  conditions,  the  toughness  of  the 
reinforced  solid  is  given  by  Rose  (1975)  as 
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Conversely,  if  /{’Cp,r/A'™a?  7,  the  results  are  closely  approximated  by  the  simple  relation 
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Equation  (4.1)  applies  for  all  ratios  of  R/L,  but  (4.2)  applies  only  for  R/L  <  0.25 


Two  important  conclusions  may  be  drawn  from  the  theoretical  results.  Firstly,  crack 
trapping  and  bridging  are  very  effective  toughening  mechanisms.  While  the  toughness 
enhancement  due  to  crack  tip  microcracking  or  crack  deflection  is  modest  at  best,  crack 
trapping  combined  with  bridging  can  improve  toughness  by  a  factor  of  5  or  more.  Secondly, 
the  p^rticlp  toughness  must  exceed  that  of  the  matrix  by  a  considerable  margin  if  pinning 
particles  are  to  be  formed.  If  <  3,  the  particles  are  always  penetrated  by  the 

crack  front.  Even  if  A'Par/AT™at  ~  15,  only  10  rows  of  obstacles  remain  intact  in  the  wake 
of  the  crack. 
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It  is  reasonable  to  assume  that  particles  remain  perfectly  bonded  to  the  matrix  in 
glasses  reinforced  by  metallic  particles,  and  the  theoretical  results  are  in  good  agreement 
with  the  experimental  data  of  Krstic  et  al  (1981)  on  reinforced  glass.  However,  the  assump¬ 
tion  of  perfect  pinning  does  not  predict  accurately  the  toughness  of  monolithic  ceramics. 

In  materials  where  the  particles  axe  weakly  bonded  to  the  matrix,  the  particles  remain 
intact  and  are  pulled  out  in  the  wake  of  the  crack.  This  has  been  modelled  by  assuming 
that  the  particles  remain  perfectly  bonded  to  the  matrix  until  the  pinning  force  exceeds  a 
critical  value  Po,  whereupon  the  particles  start  to  separate  from  the  matrix.  During  pull¬ 
out,  the  pinning  force  decreases  linearly  with  crack  opening  displacement.  Under  these 
conditions  the  toughness  enhancement  due  to  trapping  remains  the  same  as  for  perfectly 
pinned  particles,  but  the  toughening  due  to  bridging  is  reduced.  The  predicted  R-curve 
for  a  typical  ceramic  is  shown  in  Fig.  15.  The  improvement  in  toughness  due  to  frictional 
bridging  is  only  of  the  order  of  20%,  although  over  500  rows  of  grains  remain  intact  in  the 
bridged  zone.  The  R-curve  behaviour  of  the  reinforced  solid  is  strongly  dependent  on  the 
strength  of  the  bond  between  particles  and  matrix.  If  the  bond  is  strong,  the  toughness 
increases  rapidly  as  the  crack  advances,  and  the  maximum  toughness  is  reached  when  only 
a  few  bridging  particles  form  in  the  crack  wake.  In  contrast,  when  particles  are  pulled 
out  in  the  crack  wake,  the  R-curve  increases  only  slowly  with  crack  advance.  Under  these 
conditions  several  hundred  bridging  particles  may  form  before  the  maximum  toughness  is 
reached:  this  requires  a  large  crack  extension. 

The  calculated  toughening  due  to  frictional  crack  bridging  is  in  good  agreement  with 
the  experimental  measurements  in  alumina  by  Vekinis  et  al  (1990).  The  influence  of  crack 
trapping  and  bridging  can  account  for  the  discrepancy  between  surface  energy  and  fracture 
toughness  observed  in  polycrystalline  ceramics.  However,  a  key  issue  has  not  been  addressed 
by  the  present  analysis.  We  have  assumed  that  some  grains  in  the  ceramic  are  tougher  than 
others:  the  origin  of  the  tough  regions  is  not  clear.  There  are  several  possible  explanations. 
Compressive  residual  stresses  are  set  up  in  polycrystalline  ceramics  due  to  the  thermal  and 
elastic  anisotropy  of  the  grains.  Evans  et  al  (1977)  have  suggested  that  grains  in  a  state 
of  residual  compression  are  less  likely  to  be  fractured  by  a  growing  crack  and  may  be  left 
behind  as  bridging  particles.  A  second  possibility  is  that  the  crack  may  be  deflected  along 
the  interface  between  grains  and  remain  trapped  inside  the  interface.  These  are  promising 
areas  for  future  research. 
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Appendix  A 

This  section  derives  the  expressions  for  increments  in  the  influence  functions  6H(s,x ) 
and  6D(s,t)  given  in  equations  (2.32)  and  (2.33).  Rice  (1987)  gives  a  detailed  discussion 
of  the  procedure  necessary  to  update  these  functions.  Formally,  if  the  crack  front  C  is 
advanced  by  an  amount  6a(A)  to  a  new  configuration  C,  then  8H  and  8D  follow  from: 


SB(s,x)  =  ± 

f  H(A,  x)D(s,  A)<5a(A)  dA, 

Jc 

(Al) 

SD(s,t)  =  ^ 

/  D(s,  A)D(A,  t)8a(A)  dX. 

Jc 

(A2) 

However,  D(s,t )  varies  a s  (t  —  s)~ 2  a s  t  —>  s.  Consequently,  the  integrals  in  (Al)  and 
(A2)  axe  undefined  unless  <5a(A)  — ►  0  at  A  — ►  s,  in  which  case  they  are  defined  in  the 
principal  value  sense.  Rice  (1987)  has  pointed  out  that  the  integrals  may  be  regularised  for 
a  general  perturbation  by  exploiting  the  known  behaviour  of  the  crack  under  a  self-similar 
expansion,  rotation  and  translation.  The  necessary  steps  have  been  discussed  by  Bower 
and  Ortiz  (1990). 


In  order  to  apply  (Al)  to  a  general  crack  advance,  it  is  first  necessary  to  apply  a 
combination  of  a  self-similar  expansion  q,  a  translation  8x  and  a  rotation  8u>  to  the  crack 
front  so  as  to  set  up  a  reference  configuration  Cre*.  The  reference  configuration  must  be 
chosen  so  that  the  point  x  is  common  to  both  C  and  Cref ,  but  the  crack  advance  from  Cref 
to  C  satisfies  8aTef(s)  =  0.  It  is  straightforward  to  show  that  an  appropriate  combination 
of  expansion,  translation  and  rotation  are  given  by: 


„  _  ,  [x-  z(s)]  •  n(s) 

[x  -z(s)]  •  [x-z(s)]’ 

r.  .  _  -*)  x  k)-n(s) 

-  '“W  ■[*-«.)]  ■[*-*)]' 

6x  =  (1  —  q)x  —  xx  ktfu;, 


(A3) 

(A4) 

(A5) 


to  first  order  in  8a(s).  The  resulting  reference  crack  front  displacement  is 


8aTef(A)  =  8a(\)  —  [(<?  —  l)(z(A)  —  x)  4-  8u( z(A)  —  x)  x  k]  •  n(A). 


(A6) 


In  setting  up  the  reference  configuration,  i7(s,x)  changes  by  an  amount 

1 


8H(s,x)  =  H(s,x ) 


j3/2 


- 1 


2  [x  -  z(s)]  ■  [x  -  Is)] 


t 


(A7) 

(AS) 
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since  H(qs,qx )  scales  as  H(s,x)/q 3/2  under  uniform  expansion,  and  is  invariant  to  transla¬ 
tion  and  rotation.  Adding  this  contribution  to  equation  (Al),  replacing  <5a(A)  with  £are{(A) 
and  performing  some  simple  vector  manipulations  on  the  result  leads  to  the  expression 
quoted  in  equation  (2.32). 

A  similar  procedure  may  be  used  to  extend  (A2)  to  general  crack  advance.  This  time, 
it  is  necessary  to  set  up  a  reference  configuration  that  satisfies  6aTe{(s )  =  6aTe{(t)  =  0. 
Bower  and  Ortiz  (1990)  give  the  necessary  expansion,  rotation  and  translation  as 


[<5a(t)n(t)  -  £a(s)n(s)]  •  [z (t)  -  z(s)j 

9  '  '  [*(t)  -  z(s)]  •  [*(<)  -  z(s)] 

(A9) 

[i5a(s)n(s)  —  (fa(t)n(t)]  •  [(z(s)  —  z(t))  x  k] 

[z(s)  -  z(/)J  •  [z(s)  -  z(<)] 

(A10) 

Sx  =  6a(t)n(t)  —  (g  —  l)[x(t)  —  x(t)]  x  k Suj, 

(All) 

and  the  resulting  reference  crack  front  displacement  is 

<5aref(A)  =  6a( A)  -  [(, q  -  l)z(A)  +  Sx  +  z(A)  x  k<M  •  n(A).  (A12) 


Since  D(qs,qt )  scales  as  D(s,t)/q 2  under  uniform  expansion,  and  is  invariant  to  rotation 
and  translation,  the  change  in  D(s,t)  in  setting  up  the  reference  configuration  is 


6D(s ,  t )  =  D(s,  t ) 


—2  D(s,t) 


-  1 


[£a(t)n(t)  -  <5a(s)n(s)j  •  [z (t)  —  z(s)] 
\z(t)  -  z(s)j  •  [z(<)  -  z($)] 


(A  13) 
(A14) 


Including  this  contribution  in  equation  (A2),  replacing  6a  by  6aTe{  and  simplifying  the 
result  leads  to  the  expression  given  in  equation  (2.33). 
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Appendix  B 


In  this  section,  the  steps  necessary  to  derive  equations  (2.35)  to  (2.38)  from  (2.29)  to 
(2.32)  are  outlined. 

Finding  6U(x)  (equation  (2.35))  is  straightforward.  From  (2.29),  we  may  write 

6U(x)=S-^U(x)+1-  J°°  K0(s)H(s,x)6a{s)ds,  (£1) 

where  the  first  term  gives  the  contribution  to  8U  due  to  a  change  in  the  remote  load,  and 
the  second  term  gives  the  effect  of  a  change  in  the  geometry.  Since  the  crack  geometry  and 
the  crack  advance  axe  periodic,  this  may  be  re-written  as  an  integral  over  one  wavelength 
of  the  crack 

1  W2  ±2? 

6U(x)  =  -  V  K0(s)H(s  +  mc,x)6a(s)ds.  (B 2) 

2  d  —  c/2  m±?oo 

Finally,  noting  that  H(s  +  me,  x)  =  H(s,x  —  mlj)  and  using  (2.18),  we  deduce  that 


1  yc/2 

6U(x)==-  K0(s)H(s,x)6a{s)ds.  (£3) 

2  J— c/2 


The  expression  for  8Ko{s)  follows  in  much  the  samp  wav.  Begin  by  writing: 

&  x°°  r°° 

6K0(s)  =  +  /  K0(\)D(s, \)[6a(X)  -  *a(s)n(A)  •  n (s)}  d\.  (£4) 

■**  J  —  oo 

This  may  then  be  expressed  as  an  integral  over  one  wavelength  as 
c  k°°  /■+c/2  00 

6Ko(s)  =  - - f-  /  K0(\)  V  D{\  +  mc,t)[6a{\)-6a{s)n(\)-n(s)]d\.  (Bo) 

A  d-c/2  m±^oo 

Finally,  using  (2.39),  this  becomes 

8Ko(s)  =  6-^-K0(s)  +  [C/ 2  Ko(s) D(A,  s)[6a( A)  -  6a(s)n( A)  •  n(s)]  <fA,  (£6) 

A  00  Z7T  J_c/2 


Consider  next  updating  H(s,x).  Using  the  definition  (2.19),  <5H  may  be  expressed  in 


terms  of  6H : 


+  oo 

m(s,x)=  y  6H(s,x  +  mLj). 

m  =  —  oo 


(B  7) 
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Now,  SH  is  given  by  (2.32).  It  is  convenient  to  split  (2.32)  into  a  periodic  and  a  non-periodic 
part,  as 

6H(s,x)  =  [  H{\,x)D{\,s)[6a{\)  -  Sa(s) n(A)  •  n(s)] 

27T 

4-  Sa(s)'$>(s,x),(B8) 


where 


,x)  =  -L/  H(\,x+)D(\,s) 

|n(A)-n(s)- 


-  z(A)j  •  [x  -  z(s)] 


z(s)]  •  [x  -  z(s)] 
[x  -  z(A)|  x  [x  -  z(s)|  •  k 


(n(A)-n(s)] 


[x  -  z(s)]  •  [x  -  z(s)] 


[n(A)  x  n(s)]  •  k  >  d\ 


,  3„,  [x  -  z(s)]  •  n(s) 

2  '  [x-z(s)]  -[x-z(s)]' 


(B  9) 


Using  (B7),  and  rewriting  the  integrals  over  a  single  wavelength,  the  expression  for  <5H 
becomes 

,  rc/2  +  °°  +oo 

«H(S,X)  =  -  /  V  £  H( A  4-  nc,x  +  mL})D(\  +  nc ,  s)[<5a(A)  —  <5a(s)n(A)  ■  n(s)]  d\ 

~  ^  —c/2  n  =  _ go  m~  —  oo 

-f  oo 

+  <5a(s)  ^2  ^( s,x  +  m£j ).  (BIO) 

m=  —  oo 

Bv  defining  two  new  functions: 

OO 

D(s,  A)  =  £  D(\  +  nc,  s),  (B 11) 

fl—  —  OO 

-f  oo  -f-oo 

^(s,x)=  L  ^(s  4-  nc,  x  4-  mZj),  (B12) 


n=-oom=-oo 


and  writing  the  integral  in  (B9)  as  a  sum  over  successive  wavelengths,  we  obtain  the  result 
given  in  (2.38). 

Finally,  it  is  necessary  to  find  a  procedure  for  updating  G(x*n\ x**'1).  Begin  with  the 
definition  in  (2.14) 

(  1 


G(x("\x(fc))  =  { 


4R  +  TTm+  E  C(x‘">  +  mLj,x<">), 

2n2x1  m  *=  —  oo  n  =  k\ 

myt  0 


V  G(x{h)  +mL},x(n)), 


n  ^  k. 


(B 13) 


V  TO  —  —  OO 
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Each  term  in  (B13)  may  now  be  updated  in  turn.  Consider  first  the  case  n  ^  k.  From 
(2.30),  we  immediately  obtain 

'^00  i  r  +  cc 

6  G(x(n\x(k))  =  Y,  o  H(\,xik)  +mL))H(\,x{n))6a{\)d\,  n  ^  k.  (£14) 

m  =  — oo  2 

which,  written  as  an  integral  over  one  wavelength  of  the  crack  is 

1  /•+c/2 

8G{x{n\x(k))  =  -  H(A,  x^)H(A,  x(n))<5a(A)  d\,  n^k.  (B 15) 

2  J  —  c/2 

Consider  now  the  case  n  =  k.  The  infinite  sum  may  be  updated  using  a  similar  procedure 
as  for  the  case  n  ^  k,  leaving  the  term  1/4 R  +  1  / (2tt2 x[n^ ) .  To  update  this  contribution, 
it  is  necessary  to  return  to  the  original  derivation  of  G.  Recall  that  this  term  represents 
the  approximate  solution  to  the  integral  equation 


0  =  U(x)  +  f  G(t,x)p(£)d£ 
Jo 


(B 16) 


which  is  an  expression  for  the  the  pressure  required  to  pin  a  single  circular  region  on  the 
crack.  To  update  the  solution,  we  therefore  need  to  solve 

0  =  U(x)  +  8U(x)  4-  /  (G((,x)  +  fG((,x)}(p(()  +  W(M.  (BU) 
Jo 

Neglecting  second  order  terms  and  using  (B 16),  this  becomes 


0  =  6U(x)  + 


/  <5G(£,x)p(£) 
Jo 


+  G(Z,x)6p{£)dZ. 


(B  IS) 


Substituting  for  6G  from  (2.30)  then  gives 


0  =  8U(x)  + 


Jo  I  2  Jc 


H(s,£)H(s,  x)6a(s)  ds  p(£)  +  G(£.  x)<5p(0  <*£• 


In  the  spirit  of  the  analysis  presented  in  section  1,  we  now  assume  that  R  <<  L  and  write 


H(S'£)  %  H(s,x), 


accurate  to  lowest  order  in  R/L.  This  leads  to  the  following  integral  equation  for  6p(£) 

0  =  8U{x)^p\  [  H(s,£)H(s,x)6a(s)  ds  4-  /  G(£  x)6p(f)  d&  (B 21) 

3  Jc  Jo 

where  P  is  the  resultant  of  the  pressure  distribution  p(£),  as  defined  in  section  2.  Compari¬ 
son  with  (Bio)  shows  that  the  solution  to  this  equation  must  have  the  form  given  in  eqution 
(2.8).  The  distribution  of  pressure  on  each  pinning  particle  therefore  remains  constant  as 
the  crack  front  is  perturbed,  only  its  magnitude  is  affected.  Furthermore,  the  second  term 
in  (B21)  clearly  represents  the  change  in  the  1/4 R  +  \/2tt2xY  term  in  G.  Conveniently, 
combining  terms  in  6G  shows  that 

1  re/ 2 

8G(x{n\x(k))  =  -  /  H(s.x(n))H[s,x{k))8a(s)ds.  (B 22) 

2  J-c/2 
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TABLE  1 


Parameters  for  frictional  pull-out  model 


Shear  Modulus  /i 

124  GNm-2 

Poissons  Ratio  v 

0.21 

Inherent  toughness  K™*1 

1.95  MNm-3/2 

Grain  Size  2 R 

6.0  /im 

Height  of  grain  <5crit 

6/zm 

Tough  Grain  Spacing  L 

23  /im 

Residual  stress  ar 

120  MNm"2 

Friction  Coefficient  / 

0.6 

P0  =  fcTr^RScrh 

0.008  N 

A  =  (1  -  v)j[L 

6.37xl0-12  m2N-1 

P0<5crit/I2A(A'fat)2 

4.2 

R/L 

0.15 

L26CIit/ARP0 

20  x  103 

Scrit/R 

2.0 
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Fig.  1:  Mechanism  for  forming  bridging  particles  in  the  wake  of  a 
crack. 


Fig.  2:  Semi  infinite  crack  pinned  by  parallel  rows  of  obstacles, 
showing  notation  and  sign  convention. 
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Fig.  6:  The  variation  of  remote  load  with  average  crack  advance 
for  the  crack  shown  in  Fig.  4. 
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Fig.  10:  Effective  toughness  of  a  material  reinforced  by  an  array  of 
perfectly  bonded  particles.  Numbers  in  parentheses  indi¬ 
cate  number  of  stable  pinning  particles  in  the  crack  wake. 
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Fig.  12.  Variation  of  load  with  average  crack  advance  as  a  semi- 

infinite  crack  bypasses  9  rows  of  perfectly  pinned  particles. 
L/R  =  10. 
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Fig.  13:  Variation  of  stress  intensity  factors  on  bridging  particles  in 
the  wake  of  a  crack,  during  formation  of  9  rows  of  pinning 
particles.  LjR—  10 
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Fig.  14:  Comparison  of  measured  and  predicted  toughness  of  a  brit¬ 
tle  matrix  reinforced  by  tough  particles.  Theory  assumes 

^par/^mat  _  U 


Fig.  15:  Comparison  of  measured  and  predicted  R-curve  behaviour 
for  Alumina.  Parameters  used  in  the  theory  are  listed  in 
Table  1. 
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Fig.  16:  Variation  of  toughness  enhancement  due  to  frictional  bridg¬ 
ing  with  energy  dissipation.  L26CTa/ARPo  =  20  x  103, 

<5crit/P  =  2.0. 
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Fig.  17:  Variation  of  length  of  bridged  zone  with  frictional  energy 
dissipation.  L26cnt/ ARPq  =  20  x  itr,  dCTll/R  =  2.0. 
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Fig.  18:  Variation  of  length  of  bridged  zone  with  compliance  ratio. 
PoWA(A'cm“)2  =  2.0,  6crit/R  =  2.0. 


